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ABSTRACT 

A goal of forthcoming imaging surveys is to use weak gravitational lensing shear measurements to constrain 
dark energy. A challenge to this program is that redshifts to the lensed, source galaxies must be determined 
using photometric, rather than spectroscopic, information. We quantify the importance of uncalibrated pho- 
tometric redshift outliers to the dark energy goals of forthcoming imaging surveys in a manner that does not 
assume any particular photometric redshift technique or template. In so doing, we provide an approximate 
blueprint for computing the influence of specific outlier populations on dark energy constraints. We find that 
outlier populations whose photo-z distributions are tightly localized about a significantly biased redshift must 
be controlled to a per-galaxy rate of 1 - 3 x 10~ 3 to insure that systematic errors on dark energy parameters are 
rendered negligible. In the complementary limit, a subset of imaged galaxies with uncalibrated photometric red- 
shifts distributed over a broad range must be limited to fewer than a per-galaxy error rate of F cat < 2-4 x 10~ 4 . 
Additionally, we explore the relative importance of calibrating the photo-z's of a core set of relatively well- 
understood galaxies as compared to the need to identify potential catastrophic photo-z outliers. We discuss 
the degradation of the statistical constraints on dark energy parameters induced by excising source galaxies at 
high- and low -photometric redshifts, concluding that removing galaxies with photometric redshifts z ph > 2.4 
and z ph < 0.3 may mitigate damaging catastrophic redshift outliers at a relatively small (< 20%) cost in statis- 
tical error. In an appendix, we show that forecasts for the degradation in dark energy parameter constraints due 
to uncertain photometric redshifts depend sensitively on the treatment of the nonlinear matter power spectrum. 
In particular, previous work using PD96 may have overestimated the photo-z calibration requirements of future 
surveys. 

Subject headings: cosmology: theory - dark energy - galaxies: distances and redshifts - galaxies: photometry 
- gravitational lensing: weak 



1. INTRODUCTION 



Weak gravitational lensing of galaxies by large-scale struc 
ture is developing into a powerful cosmological probe (e.g. 
Hoekstra et al.|2002[|Pen et al.|2003HJaryis et al.|2003] |Van 



Waerbeke et~aTl|20051 |Jarvis et al ||20061 JS emboloni et al. 
2006 1 |Kitching et al.|2007[|Benjamin et al.|2007[|Dore et a" 
2007|]Fu et al.||2008[ l. Forthcoming imaging surveys such 
as the Dark Energy Survey (DES), the Large Synoptic Sur- 
vey Telescope (LSST), the European Space Agency's Euclid, 
and the Joint Dark Energy Mission (JDEM) expect to exploit 
measurements of weak gravitational lensing of distant source 
galaxies as one of the most effecti ve means to constrain the 
properties of the dark energy (e.g., |Hu"& Tegmark 19 99} |Hu 
Huterer|2002l peavens 2003 ; Refregier 2003 ; Refregier 



1999 



eTaLl2 004; Song & Knox 2004; Takada & Jain|20 04; Takada 
& White||2004HDodelson & Zhang||2005| |Ishak| 12005 [ |A1- 
brecht et al.|2006| |Zhan|20061 |Munshi et aL|2008j |Hoekstra 
& Jain|2008i |Zentner et al.|2008| |Zhao et al.|2009|>. The most 



stringent dark energy constraints can be achieved when source 
galaxies can be binned according to their redshifts, yielding a 
tomographic view of the lensing signal. Among the contri- 
butions to the dark energy error budget will be the error in- 
duced by the need to use approximate redshifts determined 
from photometric data (Bolzonella et al. 2000 Collister & 
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Lahav 2004 Feldmann et al. 2006 Banerji et al. 2008 ; Bram- 
mer et al. 12008: Ovaizu et al. 2008 ILima et al. 12008: Dahlen 



mer et al.|20b8[|Oyaizu et al.|2008||Lima et al.|2008[|bahlen 
et al.|20081 |Abdalla et al.|20081 |Newman||20081 |Ilbert et al. 
|2009| |Coupon et al.||2009[|Cunha et al.||2009[ |Schulz||2009| l 
because it is not possible to obtain spectroscopic redshifts for 
the large numbers of source galaxies needed to trace cosmic 
shear. Photometric redshifts are and will be calibrated by 
smaller samples of galaxies with spectroscopic redshifts. In 
this paper, we study the influence of poorly-calibrated photo- 
metric redshifts for small subsets of the galaxies within imag- 
ing samples on dark energy constraints. 

The influence of uncertain photometric redshifts (photo-z 
hereafter) on the dark energy program has been studied by 
a number of authors dMa et al. 2006; Huterer et al. 112006 



Lima& Hu 200 7~||Kitching et al.|200 8 ; Ma & Bernstein p008 
Sun et al.||2009| |Zentne~& Bhatta charya||2009| |Bernstein &| 



Huterer 2010, Zhan g et al.||2009J ] Studies of the require- 



ments for photo-z accuracy have assumed relatively simple 
forms for the relationship between the inferred photo-z of a 
galaxy and its spectroscopic redshift, in particular that this 
is a Gaussian distribution with a redshift-dependent bias and 
scatter. The underlying assumption is that this distribution 
can be calibrated with an appropriate spectroscopic sample 
over the range of redshifts of interest. These studies indicate 
that roughly N spec ~ 10 5 spectroscopic redshifts are needed 
to render photo-z uncertainty a small contributor to the dark 
energy error budget, but any particular number depends upon 
the many details of each study. Broadening the description 
of the photo-z distribution to a multi-component Gaussian 
leads to slightly more demanding requirements on the spec- 
troscopic calibration sample ( |Ma & Bernste in 2008 ). How- 
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ever, complexity or multi-modality of the photo-z distribu- 
tion will not induce large systematic errors on dark energy 
parameters, provided that this complexity is known and that 
we have some ability to calibrate complex photo-z features us- 
ing spectroscopic galaxy samples (a non-trivial assumption). 
That is not to say that dark energy constraints are insensitive 
to such complexity. Broad or multi-modal photo-z distribu- 
tions will provide an effective limit to the redshift resolution 
of tomographic weak lensing and will degrade dark energy 
constraints. If this complexity can be diagnosed in spectro- 
scopic samples, it may be treated by generalizing the model- 
ing in |Ma et al.| ( p006] > and |Ma & Bernstein] ( [2008] ). In ap- 
proximate accordance with the prevailing nomenclature, we 
refer to the galaxies for which spectroscopic calibration of the 
photo-z distribution will be possible as the core photo-z dis- 
tribution. |Ma & Bernstein| ( |2008| l studied multi -modal core 
photo-z distributions in some detail. 

These studies assume that the spectroscopic samples that 
will be obtained will suffice to calibrate the photo-z's of all 
galaxies utilized in the weak lensing analysis. However, spec- 
troscopic calibration samples may well be deficient in spec- 
tra of some subset of galaxies that otherwise may not be eas- 
ily identified and removed from the imaging sample (for ex- 
ample, see Newman| 2008 for a discussion). Consequently, 
some fraction of galaxies in forthcoming imaging samples 
may not have photo-z's that are well calibrated spectroscopi- 
cally and may have photo-z's that differ markedly from their 
true redshifts. Including such galaxies in weak lensing anal- 
yses would lead one to infer biased estimators of dark energy 
parameters. These systematic offsets in dark energy parame- 
ters may be considerable compared to statistical errors. 

We refer to such subsets of galaxies that are not well cali- 
brated by spectroscopic samples and which have photo-z dis- 
tributions that differ markedly from the photo-z distributions 
of the core galaxy samples as catastrophic photo-z outliers. 
Our chief aim in this study is to estimate the biases induced 
on dark energy estimators by catastrophic photo-z outliers for 
a variety of possible manifestations of catastrophic outliers, 
and to estimate the level at which such outliers must be con- 
trolled in order to mitigate dark energy biases. 

We consider two broad classes of catastrophic outliers, 
differentiated by the breadth of their photo-z distributions. 
We emphasize that our definition of a catastrophic outlier is 
more inclusive than previous usage (compare to |Bernstein &| 
Huterer [2010 1. Outlier populations with photo-z's that are 
confined to a small range of highly-biased redshifts make up 
the class we refer to as localized catastrophes. As an example, 
such outliers may correspond to galaxy populations in which 
spectral features have been misidentified in broadband photo- 
metric observations; the prevailing usage of the term catas- 
trophic error closely resembles our usage of the term local- 
ized catastrophe. The second class of outliers we consider, 
which we refer to as uniform catastrophes, have photomet- 
ric redshifts that are relatively unconstrained. This class may 
more naturally be associated with a level of spectroscopic in- 
completeness yielding a population of imaged galaxies with 
little information on the reliability of their photometric red- 
shifts. 

We describe our modeling techniques in § [2] We detail 
our results on the potential importance of catastrophic photo-z 
outliers in §[3] This section includes a brief discussion of mit- 
igation strategies in which we explore the possibility of elimi- 
nating subsets of galaxies in order to reduce biases at the cost 
of increased statistical errors. We discuss the implications of 



our results in §|4]and summarize our work in §|5] 

We include in this study an appendix that may help in 
comparing published results on photometric redshift calibra- 
tion requirements. All treatments of dark energy constraints 
from weak lensing rely on some approximate treatment of 
the growth of structure in the nonlinear regime. Several ap- 



proaches are in common use (|Scherrer & Bertsch inger 



Peacock & Dodds 



1996; Seljak 2000; Ma & Fry 2000; 



cimar ro et al.|2001||Cooray & Sheth 2002; Smit h et al. 



1991 



Scoc- 
2003) 



and additional parameters have been introduced to model 
baryonic processes ( |Rudd et~aL 2007 Zentner et al.||2008] 
|Guillet et al.||2009] >. In the main body of our p aper, we use 
the fitting form provided by Smith et al. (2003 1. We demon- 
strate in the appendix that estimates of photo-z calibration 
requirements depend upon the modeling of nonlinear power. 
Implementing the |Smith et al.| ( |2003[ l relation for nonlinear 
power results in significantly reduced photo-z calibration re- 



quirements as compared to previous results (e.g., |Ma et al. 
2006) that employed the Pea cock & Dodds| ( |1996| l approxi- 
mation. 

2. METHODS 

In this section we describe the methods used in our analysis. 
We begin in § |2.1| with a discussion of our treatment of photo- 
metric redshifts, including both the core photometric redshift 
distributions as well as catastrophic outliers. In § |2.2| we 
describe our weak lensing power spectrum observables. We 
describe cosmological parameter forecasting in § |2.3| and con- 
clude with a description of our fiducial cosmology and repre- 



sentative surveys in § 2.4 



2.1. Photometric Redshift Distributions of Source Galaxies 

We characterize the distribution of photometric redshifts 
through the probability of obtaining a photometric redshift 
z ph , given a galaxy with spectroscopic (or "true") redshift z, 
P(z ph \z). The distribution of true redshifts of galaxies in a 
photometric bin labeled with index i is 



high 



m(z) = n(z) dz ph P(z ph \z), 



(1) 



where n(z) is the number density of source galaxies per unit 
redshift z, n,(z) is the number density of sources per unit red- 
shift that are assigned to the z* photo-z bin, and zJ ow and z^ lsh 
delineate the boundaries of the i* tomographic bin. 
We model the overall galaxy distribution via 



n(z) oc z l exp[-(z/z ) 12 ], 



(2) 



where zo is determined by specifying the median redshift of 
the survey and the powers of redshift are representa tive of 
the distributions of observed high-redshift galaxies (Newman 
et al.p OlO). The normalization of the overall galaxy distri- 
bution is determined by the total number of galaxies per unit 
solid angle, 



N A = 



d z n(z), 



and we designate the number of galaxies per solid angle in 
any photo-z bin as 



dzni(z). 
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2.1.1. The Core Photometric Redshift Distribution 

For the purposes of our study, we consider the core galaxy 
distribution to be comprised of galaxies with a photometric 
redshift distribution that will be well characterized through 
calibration with spectroscopically-observed galaxy samples. 
Studies using existing spectroscopic galaxy samples to pre- 
dict the photo-z distributions of galaxies in future large-scale 
image surveys indicate that the core distributions may be com- 
plicated (e. g., |Jouvel et al.1[2009} |Ilbert et al.||20"09t [Coupon 
|et al.|2009] >. A common simplifying assumption in the litera- 
ture is that the photometric redshifts of galaxies in the core 
are distributed according to a Gaussian distribution with a 



redshift-dependent mean and variance (e.g., |Ma et al. [2006 
|Ma & Bernstein|2008| l, 



Pcore(z Ph \z) = 



2tt<t- 



exp 



{Z-Z" h -T bias ) 2 



(3) 



where both cr z {z) and z blas (z) are functions of true redshift, z. 
The redshift-dependent mean and variance endow this form 
with sufficient flexibility to treat a wide variety of redshift 
distributions; however, this simple model does neglect com- 
plex features that may be present in the realized photometric 
redshift distributions of future surveys. We adopt this model 
because it is a published standard against which our results 
can be compared, and because the complexity of calibrating 
the core sample of photometric redshifts is not the primary 
aim of our work. 

We compute the functions cr ; (z) and z blas (z) by linear inter- 
polation between values tabulated at 3 1 redshift points spaced 
evenly between z = and z = 3. This choice of binning al- 
lows for maximal degradation in dark energy constraints ab- 
sent prior information about the photometric redshift distribu- 
tion of source galaxies. We treat the bias and dispersion at 
each of these redshifts as free parameters in our forecasts, so 
that there are 2 x 31 = 62 free parameters describing the core 
photometric redshift distribution. For our fiducial model, we 
take o-.(z) = 0.05(1 +z) and z bias (z) = 0. 

2.1.2. Catastrophic Photometric Redshift Outliers 

Forthcoming large imaging surveys will observe a tremen- 
dous number of galaxies. It is unlikely that accurate calibra- 
tion of every class of photometric redshift distribution will 
be made, at least in part due to the limitations of obtaining 
reliable spectroscopic redshifts (e.g., |Newman|2008) and ob- 
servations of relatively rare objects. If either the uncalibrated 
objects follow the redshift distributions of the sample of cal- 
ibrated photometric redshifts, or the uncalibrated objects can 
be identified from imaging data and removed from the sam- 
ple, they will have a relatively benign impact on the dark en- 
ergy aims of these surveys. In the former case, they present 
no systematic error because they follow the redshift distribu- 
tion of the majority of galaxies, and in the latter case they can 
be removed from the imaging sample at a small cost in sta- 
tistical uncertainty. Conversely, if a sample of uncalibrated 
source galaxies that does not follow the redshift distribution 
of the calibrated sources remains in the imaging data used 
for dark energy constraints, this could represent a significant 
additional systematic error. In approximate accordance with 
established nomenclature, we refer to subsets of galaxies that 
do not follow calibrated photometric redshift distributions and 
cannot be removed from imaging data as catastrophic photo- 
metric redshift outliers. 



In practice, it is expected that catastrophic photometric red- 
shift outliers will be present at some level in forthcoming 
imaging surveys. The prevalence of multi-modal features in 
the photo-z distributions of existing calibration samples is a 
clear illustration of the difficulty of determining galaxy red- 



shifts from photometric color s ([Oyaizu et al.||2008| [Cunha 



eFaT1|2009l |Ilbert et al.||2009} |Coupon e t al. 2009j"When a 



population of galaxies responsible for a non-trivial photomet- 
ric redshift determination appears sufficiently often in spec- 
troscopic samples, its associated photo-z error can be cali- 
brated, perhaps leading to multi-modal features in the core 
distribution. However, there will inevitably be populations of 
galaxies with photo-z degeneracies that are sufficiently rare 
so as to evade spectroscopic sampling, the spectroscopic cali- 
bration of a truly representative sample will not be complete, 
and the removal of galaxies with troublesome redshifts from 
the imaging data will be imperfect. Each of these difficul- 
ties leads to a population of outlier galaxies, with distribu- 
tions not described by the core photometric redshift model, 
that contributes a systematic error to dark energy parameter 
estimators. 

To illustrate the distinction between catastrophic outliers 
and multi-modal features in the core, consider the photo-z dis- 
tribution illustrated in Figure [T] The bulk of the galaxies in 
this distribution (black diamonds) are scattered about the line 
z = z ph . This is a population of 400 galaxies drawn from the 
Gaussian distribution of Eq. Q. There are also two "islands" 
in the distribution. The appearance of these island contribu- 
tions to the photo-z distribution is quite similar, but they are 
intended to represent photo-z errors of a qualitatively differ- 
ent nature, as discussed below. One island has (z,z ph ) coor- 
dinates (0.3,3.7) and the other has (2.0,0.8). The island at 
(z,z ph ) = (0.3,3.7), consisting of black squares, is a schematic 
representation of some subset of galaxies that give a known, 
calibrated, small probability of yielding a highly-biased pho- 
tometric redshift. This is a component of a multi-modal core 
distribution and may either be calibrated with spectroscopy or 
removed from the sample. The island at (z,z ph ) = (2.0,0.8), 
consisting of the red crosses, is a schematic representation of 
a catastrophic outlier population. These are a small subset of 
galaxies with true redshifts near z ~ 2, that yield strongly bi- 
ased, but localized, photometric redshifts. Moreover, this is a 
population that is either not identified and calibrated in spec- 
troscopic samples, or is incompletely removed from imaging 
data, so that this outlier contributes a systematic error to the 
dark energy error budget. This is the type of error that is 
our focus in this manuscript. Finally, there is a population 
of galaxies that is localized near z w 1 and spread uniformly 
across z ph . These galaxies represent another extreme of catas- 
trophic photo-z errors in that the redshifts may not be strongly 
biased, but they are poorly constrained and will contribute 
systematic errors for dark energy. 

We emphasize the distinction between P(z ph \z sp ) and the 
posterior redshift distribution for an individual galaxy result- 
ing from a photometric redshift estimation algorithm, often 
denoted as p(z ph ). Often, a single redshift estimate is assigned 
to a galaxy. In this case, each point in Fig.[T]may correspond 
to the true redshift and the estimated redshift of a galaxy or 
population of galaxies. One may utilize more of the informa- 
tion in the p(z ph ), in which case the local density of points in 
Fig. [T] may correspond to regions in which the posterior has 
non-negligible support. Our aim is to outline a set of general 
impacts induced by making large, uncalibrated photometric 
redshift errors. We use P(z ph \z sp ) to quantify these effects be- 
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FIG. 1. — A toy illustration of a multi-component photometric redshift 
distribution. The aim of this figure is to provide a convenient, schematic 
representation of the photometric redshift distributions we explore. Black di- 
amonds are galaxies in the primary peak of a Gaussian core population of 
photometric redshifts specified by Eq. j3J. Black squares are galaxies in a 
secondary peak in a multi-modal core distribution. These photometric red- 
shifts are offset from the line z = z ph , but they are a known component of the 
photometric redshift distribution, and if they are represented adequately in 
spectroscopic data they can be calibrated out. Red crosses are galaxies that 
reside in a catastrophic outlier population with significantly biased, but rela- 
tively localized, photometric redshifts. In our nomenclature, this population 
is not represented in spectroscopic calibration samples and contributes a sys- 
tematic error to dark energy parameters. The red triangles represent galaxies 
that comprise a uniform catastrophic outlier population, where photometric 
redshifts are relatively unconstrained. The labels z ca t, Az ca t, and zj*. desig- 
nate the parameters of our catastrophic photometric redshift models. 



cause this allows for a very general characterization of the 
influences of photometric redshift errors (or equally, errors 
in compressing the information contained in the posteriors). 
This is sensible because P{z ph \z sp ) can be constructed from 
the posteriors of a calibration set in a straightforward manner, 
but this relationship is not invertible so that general statements 
are difficult or impossible to make. 

2.1.3. Localized Catastrophic Outliers 

One cause for localized catastrophic redshift errors (such 
as the red crosses in Fig. [TJ is the misidentification of a spec- 
tral feature in broadband photometric observations of galaxies 
over some range of true redshift. A specific example of this 
occurs when the Lyman-break is confused with the 4000 A- 
break. The effect on the photo-z distribution of a small portion 
of errors due to Lyman-4000 A confusion would look some- 
thing like the small island of squares at (z,z ph ) = (0.3,3.7) in 
Fig. [T| (Be rnstein & Huterer||20T0] >. Confusion between the 
Lyman and 4000A breaks may occur often enough in spectro- 
scopic samples to be calibrated and thus included as a sec- 
ondary peak in the core distribution, but it is possible that 
there will be other small redshift windows where tertiary is- 
lands remain uncalibrated. 

Throughout this paper we adopt a simple model for the 
photo-z distributions of localized catastrophes as Gaussians 
with spreads a cat centered away from the core at z p 



7 ph 

-cat' 



Pcat(z Ph k) = 



- exp 



(z ph - 



-7 Ph ) 2 
Zcat) 



2er*L, 



(4) 



The parameter Zc at specifies the location of the island in pho- 
tometric redshift, and er cat gives the spread of the catastrophe 
in z ph . In the presence of a localized catastrophe the total pho- 
tometric redshift distribution is 



7 ph\ 



■[l-E(z)F cat ]P core \ 

+ E(z)F cat P cat (z ph \z) . (5) 
The catastrophic error occurs over only a specified range of 
true redshifts, z cat - Az cat /2 < z < z cat + Az cat /2, as enforced 
by the function 



s(z) = e 



z-z c 



(6) 



where Q(x) is the Heaviside step function. The quantities z p at 
(location of the local catastrophe in photometric redshift), z cat 
(central value of the range of true redshifts over which the 
catastrophe occurs), Az cat (width of the range of true redshifts 
over which the catastrophic error is made), and er cat (spread 
in z ph of the catastrophe) are four of the five parameters that 
specify the local catastrophe model. The fifth parameter, F caU 
is the fraction of galaxies in the true redshift window set by 
z cat and Azcat for which the catastrophic error occurs. 

The term S(z)F C at removes the appropriate fraction 
of galaxies from the core distribution and ensures that 
J °° dz ph Ptot (z ph \z) = 1- As a concrete example, the catas- 
trophic outliers represented by the red crosses in Fig. [T] are 
galaxies drawn from our model with F cat = 0.03, z p at = 0.8, 
z cat = 2.0, cr cat = 0.1, and Az ca t = 0.1. For the sake of prag- 
matism, we present results for localized catastrophes in inter- 
esting limits of this five-dimensional parameterization rather 
than an exhaustive exploration of these parameters. 

2.1.4. Uniform Catastrophic Outliers 

Empirically, photometric redshift determination algorithms 
applied to extant calibration samples yield photometric red- 
shift estimates that are relatively unconstrainted on some sub- 
sets of galaxies. For example, the photo-z distribution of 
galaxies in both the Canada-France-Hawaii Telescope Legacy 
Survey (CFHTLS) ( |Coupon et al.|2009|) and the C osmologi- 
cal Evolution Survey (COSMOS) gilbert et al.|2009) possess 
such a feature within the range of error rates we explore in 
this work. Unconstrained photometric redshifts represent a 
regime complementary to localized catastrophes. In this case, 
photometric redshifts may be obtained with nearly equal prob- 
ability over a significant range of redshift. Such broad errors 
may occur when light from one galaxy is contaminated by 
light from another source nearby in angular separation but at 
a different redshift. It is natural to expect that such errors will 
occur most frequently near the peak of the observed galaxy 
number density n(z). 

Motivated by the presence of such errors, we also treat the 
extreme case of relatively unconstrained photometric redshifts 
by using a uniform distribution for z ph , over a symmetric win- 
dow in true redshift centered on z cat and spanning a width of 
Az ca t- We refer to this kind of error as a uniform catastrophe 
for simplicity. In the presence of a uniform catastrophe the 
total photometric redshift distribution is 

Pot {z ph \z) = [1 - H(z)F cat ] P core (z ph \z) 

+ S(z)F cat /(z 

max £min) : (7) 

where z m in and z m ax delineate the photometric redshift range 
of the survey. In analogy to localized catastrophes, the func- 
tion S(z) restricts the true redshift range over which flat catas- 
trophes occur and F ca , specifies the fraction of galaxies in this 
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true redshift window whose redshifts are catastrophically in 
error. Therefore, three parameters specify this simple model, 
namely F cat , z cat , and Az cat . The uniform catastrophe repre- 
sented by the red triangles in Fig. [T] is drawn from a model 
with F cat = 0.05, z cat = 1.0, and Az cat = 0.2. 

2.2. Cosmic Shear Tomography 

In this study, we consider constraints from weak gravita- 
tional lensing observables only. We split source galaxies into 
Atom photometric redshift bins and consider as our observ- 
ables the A^tom (A'tom + l)/2 distinct number-weighted auto- 
and cross-power spectra of convergence among the source 
redshift bins. Unless otherwise stated, we bin source galax- 
ies in equal intervals of redshift between z ph = and z ph = 3 
and take A'tom = 5, resulting in 15 distinct observables. For 
this redshift range, five-bin tomography is a useful standard 
because this binning scheme suffices to saturate dark energy 
constraints ( |Ma et al.|2006| we have verified that this remains 
so within the parameters of our study as well). 

The galaxy number count in each tomographic bin Nf , the 

cross-spectra between bins i and j, Pk(£), and the number- 
weighed spectra Vtki) are related by 



T>£(i) = NfNf !»(£)-- 



In Eq. (jSj, £ is the multipole number, H(z) is the Hubble ex- 
pansion parameter, D A (z) is the angular diameter distance to 
redshift z, and Ps(k,z) is the three-dimensional matter power 
spectrum. The lensing weight functions, W/(z), weight the 
cosmic shear signal according to the redshift distributions of 
galaxies within each tomographic bin and are defined as 

3 f°° D (z 7*) 

W i (z)=-n w Hl(l+z)D A (z) dz' * (9) 
2 J z D A (z') 

where D A (z,z') is the angular diameter distance between red- 
shifts z and z'. 

2.3. Parameter Forecasting 

We use the Fisher matrix formalism to study the constrain- 
ing power of our weak lensing observables on dark energy pa- 
rameters as well as to quantify the systematic errors on dark 
energy parameters that result from catastrophic photometric 
redshift errors. The Fisher matrix formalism is ubiquitous 
in cosmological parameter forecasting (useful references re 
lated to the present application inclu de, (Jungma n et al 



Tegmark et~aTJ[l997l |Seljak|[T997| [Kosowsky et al. 



1996 



2002 



Huterer & Takada||26051 |Albrecht et al.||20061 |Bernstein &| 
Hutere rp010| l, so we simply quote relevant results her e. The 



particular implementation we use closely mirrors that in Zent- 
|ner et al.| (2008) and |Hearin & Zentner| \2QQ9) , to which we 
refer the reader for details. 

The Fisher matrix is given by a sum over the observables. In 
the particular case of weak lensing power spectra, the spectra 
at different multipoles can be treated as independent and this 
sum can be written as 



F aP = J> + l)/ sky £ ^[C-^AB ^ + 
C,„ A,B ° Pa ° P P 



(10) 



where the "Pa are the set of observables indexed by a single la- 
bel, C -1 is the inverse covariance matrix of these observables 
at fixed multipole, [C _1 ]ab are the components of the inverse 



of the covariance matrix (we include the brackets for clarity), 
and p a are the theoretical model parameters. We choose an 
indexing scheme in which lower-case Greek letters designate 
model parameters, upper-case Latin letters designate observ- 
ables, and lower-case latin letters designate photometric red- 
shift bins, and take the mapping between observable number 
and tomographic bin number to be A = i(i- l)/2 + j. Through- 
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out this paper we use £ m i„ = 2/ sky , where / s k y 
sky coverage of the weak lensing survey, and £, 
a rough indication of the scale beyond which a number of 
weak lensing approximations break down (White & Hu 2000 
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Sembolo ni et al.|2006[|Rudd et al.|2007f 7 
The covariance matrix of observables at each multipole is 



CabW = P^(£)P£(£)+P*(£)P£(£) 



(11) 



where the indices i and j map onto A and k and I map onto B. 
The observed number-weighted power spectra, Pk(£), have 
contributions from signal and shot noise, 



PH(£) = rm+Nf6 u (ri\ 



(12) 



where the quantity (7?) is the intrinsic source galaxy shape 

noise. We conform to recent convention and fix J (7?) = 0.2, 

so that all deviations from this noise level are incorporated 
into an effective galaxy number density. 

The Fisher matrix formalism provides an estimate of the 
parameter covariance near a fiducial point in the parameter 
space. One chooses fiducial values for the model parameters 
and estimates the error on parameter a from the inverse of the 
Fisher matrix at this point, a{p a ) = LF _1 ] Q a- Within this for- 
malism, statistically-independent prior information about the 
parameters is easily incorporated by simple matrix addition. 
The second term in Eq. ( 10 1 is the prior matrix. In our analy- 
sis, we assume independent prior constraints on cosmological 
parameters, so that the prior matrix reduces to a simple diag- 
onal matrix, F*g = Sap/iCa) 2 , where 6 a p is the Kronecker-(5 
symbol and a a is the prior l-cr, Gaussian constraint on pa- 
rameter p a . We itemize our fiducial model and priors in the 
following subsection. 

Given a systematic error that induces a specific shift in the 
observables, one can use the Fisher matrix to estimate the en- 
suing systematic error in model parameters. Using APa to 
denote the difference between the fiducial observables and the 
observables perturbed by the presence of the systematic error, 
one will infer a set of parameters that is systematically offset 
from the true parameters by 



Sp a = ^F -1 W J2 (2e+ Way E APa^Jab 



(13) 

The primary results of our work are estimates of the system- 
atic errors in dark energy parameters induced by catastrophic 
photometric redshift outliers. In related literature, the 8p a are 
often referred to as biases; however, we refer to them as sys- 
tematic errors in order to avoid potential confusion with the 
biases in photometric redshifts. 

2.4. Cosmological Model and Survey Characteristics 

We assume a cosmological model specified by seven pa- 
rameters. Three of these parameters describe the dark en- 
ergy. These three parameters are the present energy density in 
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units of the critical density, f^E = 0.76, and two parameters, 
wo = -1 and w a = 0, that describe a linearly-evolv ing dark en- 
ergy equation of state, w(a) = wp+(l -a)w a (e.g.,|Linder|2003[ 
Cheval lier~& Polarski|200T] |Huterer & Turner|20011|Albrecht| 
et al.||2006| l. The values specified for these parameters are 
those in our fiducial cosmological model. In models with a 
time-varying dark energy equation of state it is interesting to 
present results for the constraint on w(a) at the scale factor at 
which it is most well constrained. The scale factor at which 
vv(fl) can be best constrained is the pivot scale factor a p , and 
is related to the Fisher matrix components as 



TABLE 1 

Representative Surveys and Baseline Constraints 



1 + 



L 1 JvVaWa 

The pivot equation of state parameter is 



and the error on vv p is 



W p = W(flp) = vv + ( 1 - fl p )w a 

([^ ]u>'oW ;1 ) 



c 2 (w p ) = [F l ] W(jW(1 



(14) 



(15) 



(16) 



The dark energy task force quantifies the constraining power 
of forthcoming surveys according to a figure of merit that re- 
flects the areas of the confidence ellipses in the wo-w a plane. 
In particular, the task forc e quotes values for the combination 
T=[a(w. d ) x cr(wp)]- 1 ( Alb recht et al.|2006 1. 

The other cosmological parameters we consider and the 
fiducial values they assume in our modeling are: the non- 
relativistic matter density cj m = il M h 2 = 0.13; the baryon 
density cjb = ^b^ 2 = 0.0223; the amplitude of the primor- 
dial curvature fluctuations A 2 ^ = 2.1 x 10~ 9 (though in prac- 
tice we vary lnA 2 ^ when computing derivatives of this pa- 
rameter) evaluated at the pivot scale k p = 0.05 Mpc" 1 ; and 
the power-law index of the spectrum of primordial density 
fluctuations n s = 0.96. We adopt relatively conservative pri- 
ors of a p (u M ) = 0.007, er p (w B ) = 10" 3 , cr p (ln A^) = 0.1, and 
cr p (n s ) = 0.04, each of which is comparable to contemporary, 
marginalized constraints on these parameters (Komat su et al.| 
2008). Using marginalized, contemporary priors allows for 
somewhat more parameter degeneracy than may be possible 
with Planck data and leads to dark energy parameter forecasts 
that are relatively conservative. 

In principle, it is relatively straightforward to scale param- 
eter forecasts from one experiment to another (e.g., Ma et al. 
|2006| [Bernstein ^fe Huterer|2010| l; however, in the interest of 
simplicity, we present explicit results for three specific ex- 
perimental configurations that span the range of observations 
expected of forthcoming instruments. 

The Dark Energy Survey is the most near-term survey that 
we consider 5 . We model a DES-like survey by assuming a 
fractional sky coverage of / s k y = 0.12 and a surface density 

of imaged galaxies of N A = 15/arcmin 2 . Second, we con- 
sider a narrow, deep imaging survey similar to a Supernova 
Acceleration Probe-like implementation of a JDEM 6,7 . We 
refer to this second type of survey as DEEP and model it with 
/sky = 0.05 and N A = 100/arcmin 2 . Lastly, motivated by a fu- 
ture ground-bas ed imaging survey such as may be carried out 
by the LSST 8 ( |Abell et al.||2009) , or a space-based mission 



http : //www . darkenerqvsurvey . org 

http : //universe . nasa . gov/program/probes/ jdem 



7 http : / /snap ■ lbl ■ gov 

8 http : //www. lsst . org 
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Survey 


/sky 


N A [arcmin 2 ] 


^■med 


a(w ) 


cr(w a ) 


CT(Wp) 


T 


DES 


0.12 


15 


0.7 


0.25 


0.77 


0.07 


18.6 


WIDE 


0.50 


30 


1.0 


0.07 


0.22 


0.02 


227.3 


DEEP 


0.05 


100 


1.0 


0.10 


0.33 


0.04 


75.6 



NOTES. — Column (1) gives the survey that motivates the particular choice of 
parameters. Column (2) is the fractional sky coverage of the survey. Column 
(3) gives the effective galaxy number density N A , in arcmin - '. We have fol- 
lowed current convention and adopted a fixed shape noise of y (7,>„) = 0.2, 
assuming deviations from this assumption to be encapsulated in the effective 
galaxy number density. Column (4) gives the median redshift of galaxies in 
the survey. Columns (5)-(8) give dark energy equation of state constraints 
in the limit of perfect knowledge of the photometric redshift distribution of 
sources. These include the uncertainty on the pivot equation of state <r(Wp) 
and the product T = [<r(w a ) X o^Wp)] -1 . Note that these constraints are from 
the weak lensing components of these surveys only and account for statistical 
errors only. 



such as the European Space Agency's Euclid 9 (Refregier et al. 
2010| l, we consider a survey with very wide sky coverage tak- 
ing / s k y = 0.5 and N A = 30/arcmin 2 . We refer to this class of 
survey as WIDE. We assume that the median galaxy redshift 
in the WIDE and DEEP surveys is z me d =1.0 and that the me- 
dian galaxy redshift in the DES-like survey is z me d = 0.7. In all 
cases, we follow recent convention by taking the shape noise 
to be \J (7?,) = 0.2, subsuming additional noise contributions 
into an effective galaxy number density. Table [1] summarizes 
our assumed survey properties. 

3. RESULTS: SYSTEMATIC ERRORS ON THE DARK ENERGY 
EQUATION OF STATE 

In this section, we present the results of our study of catas- 
trophic photometric redshift outliers. We begin with the base- 
line constraints on the dark energy equation of state param- 
eters in the limit of perfect knowledge of the source galaxy 
photometric redshift distribution in § 3.1 We continue in a 
sequence of increasing complexity. We quantify the influence 
of catastrophic photometric redshift errors in the limit of per- 
fect knowledge of the core photometric redshift distribution in 
§ |3.2| We present results on the influence of catastrophic pho- 
tometric redshift errors in the more realistic case of imperfect 
knowledge of the core distribution in § 3.3 We explore the 
prospect of excising galaxies based on their photometric red- 
shifts as a simple, first-line defense against systematic error s 
induced by catastrophic photometric redshift errors in § 3.4 



3.1. Baseline Constraints 

We begin our results section by stating our forecasts for 
dark energy constraints in the limit of perfect knowledge of 
the photometric redshift distribution. With little uncertainty 
in photometric redshift distributions, the statistical limits of 
forthcoming survey instruments would allow for constraints 
on the dark energy equation of state at the level of a few per- 
cent, as summarized in Table [T] We emphasize here that the 
limit of perfect knowledge of the photo-z distributions is not 
the assumption that photometric redshifts are precisely equal 
to the true redshifts of the source galaxies. Rather, the as- 
sumption is that there are no catastrophic errors, and that the 



http://sci.esa.int/euclid 
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photometric redshift distribution is described by the Gaussian 
in § |2.1.1 such that all 62 parameters used to specify the Gaus- 
sian distribution are known precisely. 

3.2. Systematic Errors in The Limit of Perfect Core 
Knowledge 

In this section, we present results for systematic photomet- 
ric redshift errors in the limit of perfect knowledge of the core 
distribution of photometric redshifts. This amounts to the as- 
sumption of prior knowledge of the 31 dispersion [ct z (z)] and 
31 bias [zbias(z)] parameters defined in § |2.1.1| to a level of 
< 10~ 3 , which could be achieved with a representative sam- 
ple of >4x 10 5 spectroscopic redshifts distributed in red- 
shift in a manner similar to those in the imaging survey (see 
|Ma et al.|2006| |Ma & Bernstein|2008"l and the discussion in 
the appendix of this manuscript). This is a simple case to 
begin with as it allows exploration of the influence of catas- 
trophic redshift errors over a range of the catastrophic photo-z 
parameter space without the additional complications associ- 
ated with redshift-dependent priors on the core photo-z dis- 
tribution. This is the limit explored by Bernstein & Huterer] 
d20T0l >. 

3.2.1. Uniform Catastrophes 

First, we address systematic errors induced on dark en- 
ergy parameters by a small population of uniform catastro- 
phes. Uniform catastrophes are cases in which some small 
population of galaxies with true redshifts in the range (z cat — 
Azcat/2 < z < z C at + Az cat /2 yield photometric redshift esti- 
mates that are distributed broadly in z ph . This class of error 
differs from the conventional use of the term catastrophic er- 
ror and may more naturally be interpreted as a tolerance on 
spectroscopic incompleteness. 

For simplicity, we take the central redshift of the uniform 
catastrophe to be z cat = Zmed, and determine systematic errors 
as a function of Az ca t, the width of the range of redshifts over 
which such errors occur, and F cat , the fraction of galaxies in 
this range of true redshift that correspond to this type of catas- 
trophic error. We refer the reader to Eq. (|6]l and Eq. (|7]l for 
the expressions that formally define these parameters. While 
we vary these parameters independently, they are both related 
to the total number density of sources with redshifts that are 
catastrophically in error, 



Uniform Catastrophe: Systematic Error in w 



ZcM+- 



dz' n(z'), 



(17) 



where n(z) is the overall, true redshift distribution of galax- 
ies. We should expect systematic errors to increase with both 
,F cat and Az cat because higher values of either parameter result 
in a greater total number of catastrophic errors in the outlier 
population. 

In Figure|2]we have quantified the systematic errors induced 
by uniform catastrophic errors as a function of the parame- 
ters of our simple model. The curves in Fig. [2] are contours 
of constant systematic error on dark energy parameters (for 
example, |6(wo)| for Wo) expressed in units of the statistical 
error (er(wo) for wq) at points in the Az cat -,F cat plane. For each 
of the DES, Wide, and Deep surveys, the solid curves trace 
systematic errors in dark energy that are three times the sta- 
tistical errors, while the dashed curves trace systematic errors 
that are 1/3 of the statistical error. For each of the surveys 
depicted in Fig. [2] the region of catastrophic parameter space 
that is bracketed by the solid and dashed curves labeled with 



0.010 




3a - 
1a 
o/3 - 



1.0 

4z„. 



2.0 



Uniform Catastrophe: Systematic Error in w a 



0.010 









o/3 
1o ■■■ 

3o - 






'""■°e^., t 


DEEP 























































0.0 0.5 1.0 1.5 2. 

A z cal 

FIG. 2. — Systematic errors on dark energy parameters in the case of a 
uniform photometric redshift catastrophe. The horizontal axis is the width 
of the range in true redshift over which the uniform catastrophe is realized, 
Azcat. This range in true redshift is centered at z C at = Zmed for each experiment 
(Zmed = 0.7 for DES and Zmed — 1 for DEEP and WIDE). The vertical axis 
is the catastrophic error rate per galaxy within this true redshift range, F cat . 
The solid (dashed) lines show contours of constant systematic error equal 
to three times (one third) the statistical errors on each of the dark energy 
parameters. The top panel shows contours for wq and the bottom for w a . 
Chance cancellations in the biases induced by high- and low-redshift galaxies 
cause the DES to be sensitive to catastrophic errors at similar levels to the 
WIDE survey and more sensitive than the DEEP survey. 



the corresponding survey name corresponds to outliers that 
produce systematic errors which are comparable to statistical 
errors. Systematic errors are relatively small compared to sta- 
tistical errors in the regions below the dashed curves. Each 
curve plotted in Fig. 2]has been generated with a fixed value 
of z cat = z me d. For the Wide and Deep surveys z m ed = 1, so 
when Az cat = 2 the uniform catastrophes are made over the 
true redshift range < z < 2. For DES z me d = 0.7, so once 
Az cat > 1 .4 the true redshift window over which catastrophes 
are made only increases at the high-redshift boundary. 

For each of the contours of constant systematic error in 
Fig. [2] F cat decreases with increasing Az ca t- This is simply 
because increasing the redshift range over which the catas- 
trophic errors are being made (Az cat ) leads to an increased 
total number of catastrophic errors, resulting in a decreased 
tolerance to the error rate, (F cat ). Alternatively, the total num- 
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ber of catastrophic errors in an outlier population is given by 
the integral in Eq. ( fTT] ), and the contours of constant system- 
atic error roughly trace constant values of N^ t . The contours 
flatten considerably for errors that occur over a redshift range 
Az ca t > 0.4 because there will be comparably few imaged 
sources with true redshifts near z ~ or with z > 1 .5. 

This treatment of a uniform z ph catastrophe may appear 
somewhat contrived but it gives insight into a few basic re- 
sults that are important to recognize. It is clear that the utility 
of forthcoming shear surveys to constrain dark energy is sen- 
sitive to a fractionally small population of galaxies that may 
yield poorly-determined photo-z estimates. If the error is only 
relevant to galaxies that are relatively isolated in narrow re- 
gions of true redshift, for example with Az cat < 0.1, then er- 
ror rates as high as F CSLl ~ 1% in this region of true redshift are 
tolerable. This is simply because errors that occur with a fixed 
rate over a small redshift range result in a small total number 
of catastrophic outliers to corrupt the weak lensing tomogra- 
phy. On the contrary, if such an error occurs for a subset of 
galaxies with true redshifts in an interval of width Az cat > 0.1, 
then the error rate per galaxy must be significantly lower than 
^cat ^ 0.01 in order to render the systematic errors on dark 
energy equation of state parameters small. 

The limit of Az cat ^0.1 is interesting to consider. This 
may correspond to the case of a small fraction of galaxies 
that yield very poorly-constrained photometric redshifts over 
a broad range of true redshifts and that otherwise cannot be 
identified and removed from the imaging survey. In this case, 
the systematic error from catastrophic photometric redshifts 
becomes a considerable portion of the dark energy error bud- 
get at a rate of only F c . dt ~ 10~ 3 . Reducing the systematic 
error due to such an outlier population to a negligible level re- 
quires reducing the occurrence of such an outlier population 
to F cat <4x 10~ 4 . Strictly speaking, Fig. ^corresponds to 
errors that occur when the true redshift bana over which the 
uniform catastrophe occurs is centered on z ca t = z m ed, but for 
the Az cat > 1 limit, similar results hold for a wide range of 
z cat near unity, so this result is of some general relevance to 
photometric redshift calibration studies. 

We conclude this section with a discussion of cancellations 
that may occur among systematic errors. It may seem some- 
what surprising DES exhibits comparable sensitivity to uni- 
form errors as WIDE and is more sensitive than DEEP as 
shown in Fig. 12] In § |3.2.2| we will discuss systematic er- 
rors from local catastrophes. In particular, we will show that 
large biases occur for low z ca t and for higher z ca t just over the 
median redshifts of the surveys (see Fig. [3]). These biases have 
opposite signs and partially cancel in our forecasts for both the 
DEEP and WIDE surveys. DES is less sensitive to biases from 
galaxies in the low-redshift range 0.4 < z cat < 0.6 that get mis- 
placed to higher redshifts because these shifts must compete 
with the larger shot-noise of DES. The degree of cancellation 
depends upon modeling choices, such as fiducial model and 
cosmological parameters, but the occurrence of this cancella- 
tion is robust. 

3.2.2. Localized Catastrophes: Details 

Localized catastrophes correspond to the case where a 
small fraction of galaxies near some true redshift z C at yield 
photometric redshifts that are narrowly distributed about a bi- 
ased value Zc at tn at is very different from the true redshift z cat . 
Such errors could arise due to incomplete calibration by spec- 
troscopic surveys or from difficulty in removing troublesome 



galaxies from the imaged galaxy sample. A known example 
of such an error occurs when photo-z algorithms confuse the 
4000A break with the Lyman break, but other isolated islands 
of biased z ph persist in contemporary photo-z algorithms (see, 
e.g. [Coupon et al. 2009 , Il bert et al.|2009[ ) and may be relevant 
to forthcoming imaging surveys. 

The class of localized photo-z catastrophes is more com- 
plex than the uniform case because there are more relevant 
parameters needed to specify the manner in which a local- 
ized outlier population is distributed in z ph . Our toy model 



requires five parameters (see § 2.1.2 and Fig. [T] for an illus- 
tration). Two are the central value of the true redshift over 
which this error is operative (z ca t) an d the width of the true 
redshift range over which this error is operative (Az ca t)- Like 
in the uniform case, some fraction F cat of galaxies with true 
redshifts in the interval z C at _ Az ca t/2 < z < z ca t + Az cat /2 are 
catastrophically in error. The final two parameters specify the 
biased distribution of photometric redshifts that these galaxies 
are assigned. These are the (systematically erroneous) value 
of the photometric redshift Zc at , and the dispersion in the catas- 
trophic photometric redshift distribution (T cat about z^- 

We make an effort to remain agnostic about the classes of 
photo-z errors that may be realized in future imaging data. 
However, a complete mapping of even the simple parameter 
space we have specified for catastrophic photo-z's would re- 
quire a lengthy discussion, so we explore useful limits of the 
model parameters in order to distill our results into a small 
number of points. We are particularly interested in the limit 
where the source galaxies are placed in a narrow range of bi- 
ased photo-z (<7 C at "C 0.3 or so) because the limit of large dis- 
persion in the catastrophic photometric redshift population is 
similar to the uniform catastrophe of the previous section. 

We first isolate the sinister regions in the space of ZcarZcat 
that lead to the most destructive systematic errors in dark en- 
ergy parameters. At a set of points in the parameter space of 

(Zcat,Zcat)> we nave calculated the systematic error induced in 
vv'o and vv a by distributing some fraction F C!it of the galaxies 
with true redshifts near z ca t in photometric redshifts centered 
around some Zc at that is generally very different from z cat . We 
sample a range of values of true redshifts from z ca t = 0.05 
to Zcat = 2.95, evenly spaced in redshift intervals of 5z = 0.1 
and likewise for the photometric redshifts, z^. In the in- 
terest of simplicity, we fix the remaining parameters of our 
catastrophic photo-z model to F cat = 0.05, Az cat = 0.05, and 
cr ca t = 0.01 to isolate the dependence of the parameter bias 
upon the location of the catastrophe. 

It is important to note explicitly that we present results here 
at a fixed error fraction F CSLt , and a fixed true redshift win- 
dow width and Az cat . However, even with these parameters 
fixed the absolute number of errors varies with z ca t accord- 
ing to Eq. (17 1, which is roughly A^ at ~ «(Zcat)Az cat fcat for 



sufficiently small Az ca t, along the lines of the analogous dis- 
cussion for the uniform catastrophe in § 3.2. 1 The aim of this 



calculation is to map out the relative importance of making 
errors at a fixed rate per galaxy as a function of the true and 
photometric redshifts of the outliers. 

The results of this exercise are depicted in Figure|3] In each 
column of Fig. [3] there are two panels, corresponding to the 
systematic errors in wq and w a , for each representative experi- 
ment. The horizontal axes show values of z ca t and the vertical 
axes show values of z^- The systematic error is represented 
on the grid of (z C at,Zcat) by the color in each of the cells. In dis- 
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FIG. 3. — The severity of localized catastrophic errors as a function of the values of source Zcat, and target z^ t , of the catastrophic errors. Along the horizontal 

axes are the values of z cal while the vertical axes show zj^,, just as in Figure [l] Each point on this grid corresponds to a localized catastrophe with a fixed per- 
galaxy error rate of F ca t = 0.05, and fixed values of both the photo-z spread <7 cat = 0.01 as well as the width of the true redshift range over which the catastrophic 
error is made, Az cat = 0.05. The effect of these catastrophes on wg is shown in the top row of panels, while the systematic error on w a is shown in the bottom row. 
The absolute value of the induced systematic error is color coded; the numerical values labeling the color table to the right of each panel indicate the systematic 
error in units of the statistical uncertainty in the limit of perfect core calibration. 



cussing the results of this exercise, we find the terminology of 
|Bernstein & Hutererj ( |2010[ ) to be a useful, descriptive short- 
hand. We will refer to the tomographic bin that contains the 
z cat value of an outlier as the Source Bin of that catastrophic 
photo-z population. We call the bin containing z^ its Target 
Bin. This is because galaxies with true redshifts near Zcat ^ue 
erroneously placed in the Target Bin containing the redshift 
zj?a t . Our sampling guarantees that no localized outlier strad- 
dles a tomographic bin boundary so there are always unique 
Source and Target bins. Outlier populations that straddle a 
boundary dividing two tomographic bins can be substantially 
more severe than those that do not because such an outlier 
simultaneously contaminates multiple Target bins. We have 
chosen to ignore such outlier populations for simplicity, but 
such outliers can be modeled by two catastrophic outlier pop- 
ulations, one for each affected target bin. We will return to 
the issue of tomographic binning and straddling outlier popu- 
lations below. 

The prominent block-like features in Fig.|3]reflect the tomo- 
graphic redshift bins used in our analysis. The tomographic 
bins of the source and target galaxies largely determine both 
the magnitude and sign of the induced systematic error in dark 
energy parameters. This gives rise to features that reflect the 

structure of the photometric redshift binning in the (z C at,Zcat) 
plane. Indeed, for fixed Target and Source Bins, the specific 
value of the target redshift, z^,, within the target photomet- 
ric redshift bin has little influence on the severity of the sys- 
tematic error. However, small steps in z^at can lead to large 
changes in systematic error when the boundary dividing two 
tomographic bins is crossed. 

Varying the location in true redshift, z cat , leads to some- 
what more significant changes in dark energy systematic er- 



ror. Changing z ca t within fixed Source and Target Bins can 
result in up to a factor of two difference in systematic errors. 
Two factors primarily determine the severity of the systematic 
error as a function of the true redshift of the galaxies, z C at- 

The 

primary factor stems from the fact that a fixed fractional error 
rate (F c . dt ) corresponds to a different absolute number of errors 



i\Q t as a function of redshift, z ca t- This is reflected in Eq. ( 17 1 



The number of errors N^ t will be relatively large in a region 
near the median redshift of the survey, where the number of 
source galaxies per unit redshift, «(z), is largest. There are 
relatively few galaxies at low and high redshift, so for a fixed 
error rate, outlier populations with low or high true redshifts 
contribute a relatively small absolute number of galaxies with 
highly-biased redshifts. 

Secondly, an outlier population naturally results in a more 
severe systematic error the more the photometric redshift is 
biased away from the true galaxy redshift. Consider the region 
of catastrophic parameter space near (z ca t,2cat) = (1-5,2.7) in 
either color plot for DES. Catastrophes in this region of pa- 
rameter space correspond to outlier populations whose Source 
bin is the third tomographic bin and Target bin is the fifth 
tomographic bin. Outliers in this region of parameter space 
are assigned photo-z's that are significantly too high. De- 
creasing the value of z cat (the "source" of the error) increases 
the distance between core and outlier populations, thereby in- 
creasing the systematic error. This behavior contributes sig- 
nificantly to the systematic error gradient near (z ca t,Zcat) = 
(1.5,2.7) for DES in Fig. |3] 

Each of the three representative experiments that we con- 
sider has two distinct "hot spots" in Fig. [3] that correspond 
to the most severe types of error given a fixed error rate per 
galaxy, F cat . A common feature of all these hot spots is their 
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z cat location. Each of the hot spots lies at a z cat slightly be- 
yond the median survey redshift. This is sensible because for 
a fixed error rate, the absolute number of catastrophic errors is 
greatest when they are made at the peak in the overall galaxy 
distribution, that is near z cat = Zmed- The most damaging sys- 
tematic errors occur when the galaxies are shifted to either 
very low or very high photometric redshifts, when the target 
redshift, Zc at , is very different from the source redshift, z cat , 
because the galaxies in error are then placed at distances sig- 
nificantly different from their true redshifts. For our WIDE 
and DEEP surveys the largest systematic errors tend to occur 
for galaxies shifted from a source redshift z ca t near z me d = 1 to 
very low photometric redshifts. 

The pattern of the DES catastrophic photo-z "hot spots" dif- 
fers from that of the WIDE or DEEP surveys. Outliers with 
large values of z£ at , that is those with a target in the fourth 
or fifth tomographic bin, are relatively more severe for DES. 
This is driven by the (assumed) comparably low redshift ex- 
tent of imaged sources for a DES-like survey (with median 
redshift z m ed = 0.7) This renders a contamination that extends 
to high redshift more disruptive due to the small population 
of galaxies with truly high redshifts. Though less striking, it 
is also evident in Fig. [3] that the DEEP survey is somewhat 
more sensitive than the WIDE survey to contamination of its 
fourth and fifth tomographic bins. The differences here are 
likewise driven by different survey depths and sky coverages. 
A deeper, but narrower survey (a JDEM perhaps) is relatively 
more sensitive to small-scale fluctuations induced by structure 
at high-redshift, so disruptions to the higher tomographic bins 
are more statistically significant for DEEP than for WIDE. 

Finally, we return to the issue of tomographic binning with 
respect to the systematic errors in Fig. [3] We noted above 
that systematic errors in dark energy parameters can become 
markedly worse when the biased photometric redshifts (zj^) 
distribute galaxies across the boundary of a photometric red- 
shift bin. The reason is because two sets of observables, 
namely the auto and cross spectra associated with the two tar- 
get photometric redshift bins, become corrupted by the catas- 
trophic photometric redshift error. The implication is that the 
level of systematic error induced by a localized catastrophic 
error is quite sensitive to photometric redshift binning. This is 
contrary to the statistical errors, which are insensitive to bin- 
ning more finely than A^tom ~ 5 over the range < z < 3 ( |Ma| 
|et al.|2006> . 

Indeed this is the case. The general pattern shown in Fig. [3] 
is physically quite sensible and is robust to binning. How- 
ever, in the case of localized catastrophes, binning more finely 
may reduce the absolute amplitude of systematic errors if the 
catastrophes do not occur near the edge of a photometric red- 
shift bin. This is because smaller tomographic bins result in 
a smaller fraction of source galaxies that belong to a contam- 
inated bin. This may be useful because even in the absence 
of significant prior indications of a localized catastrophe, re- 
analyzing the data with different photometric redshift binning 
schemes may reveal potential local catastrophes. In the least, 
it should be a useful strategy to choose photometric redshift 
bins such that suspect regions of z ph , where localized catas- 
trophes may be anticipated, are contained in individual bins. 

3.2.3. Localized Catastrophes: Summary 

A succinct distillation of the dominant effects that deter- 
mine the structure of Fig.[3]is as follows. The systematic error 
induced by a localized catastrophe will be most severe when: 



1- Zcat ~ Zmed, which maximizes the total number of out- 
liers; 

2. the distance between z cat and z[? at is significant; 

3. and when z^t is such that the photometric redshift bins 
contain a fractionally large contaminant (in practice, 
high and low redshift extremes). 

The details governing the magnitude of systematic errors gen- 
erated by different regions of catastrophic error parameter 
space can be complicated. In general, these details depend 
on the relative statistical weights of the affected redshift bins, 
as well as the characteristics of the survey. 

In isolation, Fig. [3] is useful in identifying the redshift er- 
rors that most seriously compromise dark energy constraints. 
A shortcoming of Fig.[3]is that we have assumed catastrophic 
errors that occur at a fixed rate of F cat = 0.05 and are active 
only over a range Az ca t = 0.05. The systematic errors induced 
on cosmological parameters scale approximately with the to- 
tal number of catastrophes, N^ t , in Eq. ( 17 I. In practice, scal- 
ing the systematic errors to new values ofAzcat can be enacted 
over an interesting range of the parameter space by approxi- 
mating NL ~ F clit n(Zcat) AZcat- 

Figure [4] demonstrates the validity of scaling systematic 
error by the total number of errors, NjL, for three example 
localized catastrophes. Together, Fig. [3] and Fig. |4] provide 
a blueprint for estimating the systematic error induced by a 
wide range of localized catastrophes. One first reads off the 
systematic error level from Fig. [3] for the grid point of inter- 
est. For definiteness, suppose this systematic error in either 
of wo or w a is 5. Provided that Az ca t is small, one can ap- 
proximate the systematic error induced by a different effective 
value of Azfyt or F cat (call it 8') by scaling 5 in proportion to 



< t [Eq.(17i], 



S'^Sx (00.05) x(A4 t /0.05) 



(18) 



In § 3.2 we presented results on the influence of catas- 



trophic, uncalibrated photometric redshift errors on the sys- 
tematic error budget for dark energy parameters wq and w a . In 
that section, we assumed that the bulk of photometric redshifts 
had been well characterized by spectroscopy. In the nomen- 
clature of this and other papers, we assumed the limit in which 
the core of the photometric redshift distribution is calibrated 
so that its uncertainty does not contribute to the dark energy 
error budget. We developed guidance on how to optimally fo- 
cus photo-z calibration efforts and identified the most severe 
types of catastrophes. In this section, we drop the assump- 
tion of arbitrarily precise calibration of the core populations 
of photometric redshifts. Our goal is to assess the relative 
importance of calibrating the core photometric redshift distri- 
bution compared to eliminating catastrophic errors. 

We assume that the core photometric redshift distribution 
is specified by a Gaussian with redshift-dependent mean and 
dispersion. Following |Ma et aL| p006 ), we specify the un- 
known mean and dispersion at 3 1 points spaced evenly in red- 
shift from z = to z = 3 and allow for uncertainty in these 
parameters. In the interest of simplicity, we consider a one- 
parameter family for the prior knowledge about the core pho- 
tometric redshifts that may be provided by a spectroscopic 
calibration sample. We do this by assuming a representative 
population of iV spec galaxies with spectroscopic redshifts, dis- 
tributed evenly in redshift from z = to z = 3, which can be 
used to calibrate the core photometric redshift distribution. 
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FIG. 4. — The scaling of systematic error in dark energy parameters with the fractional number density of sources whose photometric redshifts are catastroph- 
ically in error for our WIDE survey. On the vertical axes are the absolute value of the systematic error in Wg (left panel) and w a (right panel) in units of statistical 
uncertainty. On the horizontal axis is N^ t /N A S3 Fcat»(Zcat)Azcat, where «(z c!1 t) is the overall redshift distribution of sources. For the dashed curves, Az ca t = 0.1 
and we increase N^ t by increasing F cM . These curves are all linear, as they should be. For the dotted curves F CM = 0.03 and we increase N^ t by increasing 

Az cal . These curves grow approximately as the linear, dashed curves. Three different catastrophic error localizations are color coded as (z c at,z^ t ) = (0.9,0.3) in 
red, (0.9, 2.7) in blue, and ( 1 .5 , 2.7) in green. The difference in intrinsic severity between these outlier populations is reflected by the slope of the corresponding 
curves, with the steeper lines corresponding to the more severe systematic errors. The agreement between dashed and dotted lines for each outliers demonstrates 
that the systematic errors induced by sufficiently well-localized catastrophes (Az cal < 0.3) scale approximately linearly with Az cat over an interesting range. 



We implement core calibration by introducing priors on the 
values of the dispersion and bias at the /th point in redshift. 
These priors are 



Az bias = 



1 



2N! 



spec 



(19) 
(20) 



NL 



where z' hias is the bias at the /th point in the tabulated core 
distribution, a\ is the dispersion at this redshift, and A^ pec is 
the number of spectroscopic galaxies in each of the 31 bins of 
width Sz = 0.1 used to calibrate the core photo-z redshift dis- 
tribution. This prior model is certainly simplistic. For exam- 
ple, in our analysis we have chosen for the sake of simplicity 
to set all of the iVL ec equal to each other, so that our implemen- 
tation assumes that calibrating spectra are sampled equally in 
redshift, whereas in reality we will have much looser con- 
straints on sources at high redshift than those at low redshift. 
Moreover, both core calibration and the ability to identify 
catastrophic outliers improve with larger spectroscopic sam- 
ples. However, we consider these issues independently in the 
interest of completeness because the details of how a realistic 
calibration program may proceed remain uncertain. 

3.3. Catastrophic Redshift Errors with Core Uncertainty 

Figure [5] is a contour plot depicting the systematic errors 
in wo and w a induced by the worst-case-scenario catastrophes 
determined in § 3.2 The prior core knowledge is specified 



by A^ S pec, which runs along the horizontal axis. The error rate, 
F cat , runs along the vertical axis. For uniform catastrophes the 
worst case outliers span the true redshift range of the survey. 
For localized catastrophes the most sinister outliers lie at the 
points of maximum systematic error in Fig. [3] The dashed 
(solid) curves are lines of constant systematic error at a level 
of one-third (three times) the statistical error on each param- 
eter. Clearly then, systematic errors are dominant above the 



solid curves and become unimportant well below the dashed 
curves. In the bottom panels we have included dotted curves 
to emphasize the region of parameter space where systematic 
errors are equal to statistical errors. These la contours are 
omitted in the upper panels to avoid clutter, but the linear de- 
pendence of the induced systematic error on the catastrophic 
error rate ensures that the la can be estimated by scaling the 
3cr or a/3 contours by a factor of three. 

Several aspects of Fig. [5] are worthy of note. The contours 
all become very flat at large N spec . This is the limit in which 
the core photo-z distribution is calibrated sufficiently well that 
it no longer con tributes to the erro r budget of wq and w a (e.g., 
Ma et al.||2006| |Ma & Bernstein||2008[). This corresponds to 



the limit of perfect knowledge of the core photo-z distribu- 
tion, and accordingly, the systematic errors asymptote to those 
quoted in §|3.2|at large N s?ec . 

For a fixed level of systematic error, experiments gener- 
ally become less tolerant of catastrophic outliers as N svec in- 
creases. This behavior is reflected in the negative slope at 
the low A^ pec -end of the contours of constant systematic error 
in Fig. [5] This is an explicit manifestation of the competi- 
tion between calibration of the "core" population of photo- 
metric redshifts and the ability to diagnose and eliminate a 
sub-dominant, poorly-understood "catastrophic" outlier pop- 
ulation. The reason for this is simply that systematic errors 
must be better understood for samples with smaller statisti- 
cal uncertainty. If the statistical errors in the measurement 
are intrinsically large, as they would be in the limit of poorly- 
calibrated photo-z's for the majority of the imaging sample, 
then high-rates of catastrophic outliers are tolerable because 
the systematic they contribute is not large compared to the 
statistical error induced by a poorly constrained core distribu- 
tion. 

When the core distribution is very well calibrated, most ob- 
vious at A^spec > 10 5 for uniform catastrophes in DES, the con- 
tours of constant systematic error transition to slightly pos- 
itive slope. This occurs when the core distribution has been 
sufficiently well calibrated that degeneracies between the pho- 
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FIG. 5. — Contours of constant wq and w a bias from the worst case catastrophe in units of the statistical uncertainty of the survey. Systematic errors in wq 
appear in the left panels and w a in the right panels. Results for the worst case uniform catastrophe appear in the upper panels, and were generated with Az ca i = 1.5 
and Zcat = 1.5. Contours of systematic error produced by localized catastrophes appear in the bottom panels. Each of the localized contours have been calculated 

with AZcat = 0.1 and o cat = 0.03. For DEEP and WIDE z«„ = 1.15. and zgj = 0.15, and for DES Zcai = 0.85, and z^ai =0.15, in accordance with the results 
illustrated in Fig. [5] 



tometric redshift parameters of the core distribution and cos- 
mological parameters are no longer significant. Calibrating 
beyond the level required to break degeneracies between cos- 
mology and the core photo-z parameters results in a slight re- 
duction in systematic errors on cosmological parameters. This 
is a specific manifestation of the general result that improving 
priors can only lead to a net reduction in the systematic er- 
rors of inferred parameters, a result discussed in considerable 
detail in |Bernstein & Huterer ( 2010| ). Clearly, the reduction 
in systematic error at very large N spec is not significant in the 
cases of interest here. 

To illustrate the competition between core calibration and 
the removal of outliers, consider some explicit examples. In 
the case of the uniform catastrophe, our Deep (Wide) sur- 
vey can tolerate catastrophic errors at a rate F c . dt > 1% if 
the core calibration is worse than the statistical equivalent of 
Nspec <3x 10 4 (A^spec < 6 x 10 4 ). For both surveys, even 
the worst-case, localized catastrophes can occur at a rate of 
F m > 1% if A^ S pec S 10 4 . Of course, the worst-case localized 
systematic errors are more subtle to interpret, as we have as- 
sumed they are only actively affecting galaxies over a range of 
true redshifts with width Az cat = 0.1; however, the magnitude 
of the induced systematic errors produced by localized catas- 
trophes active over different redshift ranges scales in propor- 
tion to A^ at ~ "(Zcat)Az cat F cat , as illustrated in Fig.ji] Detailed 
results are complex, but two simple conclusions are clear: 



1 . Limiting uniform catastrophic error rates to F cat <4x 
10" 4 (Feat < 2 x KT 4 ) for DES and DEEP (WIDE) will 
render them unimportant. 

2. Limiting individual localized catastrophic error rates to 
,F cat (Az cat /0. 1) < 10~ 3 will render them unimportant for 
each experiment. 

In practice, some amount of uncertainty in the calibration 
of the core distribution is inevitable, so error rates higher by a 
factor of a few may be tolerable, but in detail this will depend 
upon the nature of the error and the properties of the core sam- 
ple of well-calibrated photometric redshifts. Fig. [3] and Fig. [5] 
contain the information necessary to diagnose the systematic 
error for a variety of idealized, but interesting cases. 

3.4. Mitigating Systematic Errors by Sacrificing Statistics 

In § |3.2| and § |3.3| we estimated the systematic errors that 
could be induced by two broad families of catastrophic pho- 
tometric redshift error, remaining relatively agnostic about 
the source of the error. We found generally that error rates 
must be kept to levels below F cat ~ 10~ 3 , or one of a thou- 
sand imaged galaxies with large, uncalibrated redshift errors 
in order for systematic errors not to contribute to the dark 
energy error budget (though specific tolerances depend upon 
several details). This will be a relatively challenging goal for 
a photometric redshift calibration program to attain. DES, 
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JDEM, EUCLID and LSST will all require calibration of very 
faint galaxies, where precise photo-z's are difficult to obtain. 
Moreover, the types of galaxies imaged, and for which spectra 
may be available, varies as a function of redshift, so some un- 
derstanding of the details of galaxy evolution will be needed 
in order to achieve calibration goals. 

It is natural to explore simple methods to sacrifice some 
of the statistical power of imaging surveys in order to miti- 
gate larger systematic errors. One of the simplest techniques 
we can employ to limit the effect of catastrophic outliers is 
to place cuts on the range of photometric redshifts utilized 
to infer cosmological parameters (Bern stein & Huterer|2010| 
have explored such cuts for a particular model of photo-z out- 
liers). 10 The most damaging catastrophic errors are those that 
take galaxies near the median redshift of the survey and scatter 
them to significantly lower or higher redshifts, so it is sensible 
to explore the losses in statistical power incurred by excising 
galaxies at the low- and high-redshift ends of surveys. 

We demonstrate the utility of photometric redshift excision 
in this section by exploring a class of simple excision algo- 
rithms. In particular, we cut out all galaxies with photometric 
redshifts greater than some value, z™ax> an d smaller than some 
value z^ n - Figure [6] shows the statistical errors on wq and 
w a as a function of z^ ax and z^ n for our Wide survey, whose 
characteristics are similar to those expected from an LSST- 
or Euclid-like survey. The relative costs depend mildly upon 
survey parameters. 

Excising galaxies with photometric redshifts lower than 
z ph ~ 0.3 results in only a ~ 7% increase in the statistical 
errors on dark energy parameters. Likewise, excising galax- 
ies with z ph > 2.4 results in only a ~ 10% degradation in wq 
and w a constraints. Excising both of these regions of photo- 
metric redshift leads to a reduction in constraining power of 
<20%. Fig.|6] is a valuable itemization of the statistical losses 
incurred by redshift cuts and indicates that excising low- and 
high-redshift portions of the imaging surveys may be an ef- 
fective method to mitigate the influence of catastrophic pho- 
tometric redshift errors at little cost in statistical error. 

While Fig. [6] quantifies the cost of excising regions of pho- 
tometric redshift, the parametric complexity of catastrophic 
photo-z errors makes specific statements about the benefit of 
such cuts more difficult. In the case of a localized catastrophe 
that places galaxies erroneously in the excised high- or low- 
redshift ends of the survey, the induced bias can be nearly 
completely removed at the cost of the statistical degradation 
in Fig. [6] We have begun a preliminary study of the benefits of 
redshift excision, including the case of uniform catastrophes. 
In the case of our WIDE survey, excision can considerably 
reduce systematic errors induced even by the uniform catas- 
trophe when the core is not well-calibrate (N svec < 10 5 ), but 
this strategy is only of marginal value in the limit of a well- 
calibrated core. We limit the present discussion to the item- 
ization in Fig.[6]and relegate further study of redshift cuts and 
possible self-calibration of specific types of catastrophic error 
to a follow-up study. 

4. CONCLUSIONS AND DISCUSSION 

We have studied the potential systematic errors that may be 
induced in dark energy parameters inferred from forthcom- 
ing weak lensing surveys as a result of a population of source 



galaxies with photometric redshifts that deviate significantly 
from their true redshifts. We used a particular operational 
definition of catastrophic photo-z errors that is subtly distinct 
from the use of this term in some of the existing literature. 
Throughout this work, the term catastrophic photometric red- 
shift error refers to cases in which photo-z estimates differ 
significantly from true redshifts, the nature of the error has 
not been identified or calibrated with an accompanying spec- 
troscopic data set, and the outlier population has not been re- 
moved reliably from the imaging data prior to the construction 
of shear correlation statistics. One way to interpret our results 
is as requirements for spectroscopic calibration of outliers and 
the completeness with which outlier galaxies must be culled 
from the data set in order to render systematic errors in dark 
energy parameters small. 

In order to provide relatively general guidelines on the fi- 
delity with which outlier photo-z's must be understood, we 
have taken an agnostic position on the nature of what types 
of catastrophic photometric redshift outliers may be realized 
in forthcoming imaging data. This eliminates the need to an- 
ticipate what types of photo-z errors may occur at very small 
fractional rates in order to assess their general influence on 
dark energy parameters. To be sure, there are reasonable 
guesses that can be made regarding the nature of photomet- 
ric redshift errors and many algorithms exist that estimate 
redshifts from photometric data and refine estimates based 
upon comparisons with large, spectroscopic data sets (e.g., 



10 |Nishizawa et al.| {2010} also study the ability to employ photometric 
redshift cuts to mitigate the effects of catastrophic outliers, which became 
available on the Arxiv while we were submitting this manuscript for publica- 
tion. 



I Bolzonella et al.|2000||Collister & Lahav|200 4; Oyaiz u et al. 
2008 1 Feldmann et al. 2006 1 Brammer et al. 2008 ; Margoniner 
& Wittman|2008[|Cunha et al.|2009||Ilbert et al.|2009||Coupon 
|et al.||2009p . However, we have not adopted any particular 
template for photometric redshift outliers. Instead, we have 
studied two extreme limiting cases of catastrophic photomet- 
ric redshift error. 

In the first class of photometric redshift error, which we 
dubbed the uniform catastrophe, photometric redshifts are 
poorly constrained and scattered over a broad range (see, 
e.g., Ilbert et al. 2009; Coup orTet al.| 2009 for examples of 
such features). Photo-z errors resembling our uniform type 
must be well controlled. If such errors occur even for a rel- 
atively small fraction of galaxies near the median redshift of 
a given survey, the systematic errors induced on dark energy 
parameter estimators will be significant. Roughly speaking, 
we find that the error rate per galaxy must be maintained at 
^cat ;$ a few x 10~ 4 . However, the uniform catastrophic error 
is a relatively simple variety so that self-calibration may well 
be feasible. One could resign oneself to the fact that such an 
error will occur and add the error rate F cat (and perhaps other 
parameters such as Az cat ) to the set of nuisance parameters 
to be marginalized over. This self-calibration could eliminate 
the systematic error, but will broaden statistical errors. We ex- 
plore self-calibration of particular catastrophic photo-z errors 
in a forthcoming paper. 

The second class of errors, which we refer to as localized 
catastrophes, takes source galaxies with particular true red- 
shifts and assigns them photometric redshifts with a large bias 
but small scatter. Localized catastrophes have a broader range 
of possibilities and are more difficult to deal with. Fig. [3] 
Fig. [5] and Eq. Q constitute a blueprint for estimating the 
severity of a broad range of possible localized photometric 
redshift catastrophes. Quite generally, we find that the sys- 
tematic errors they induce are sensitive to the scheme used 
to bin the source galaxies in photometric redshift. This sug- 
gests that an iterative scheme of re-binning may be an effec- 
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FIG. 6. — The statistical cost of excising low- and high-redshift shear information on constraints of wq (left panels) and w a (right panels) for our WIDE survey. 
In the top row, the value of the maximum photometric redshift of the survey appears along the horizontal axis while the different lines show different choices of 
the minimum photometric redshift as indicated. Along the vertical axis is the fractional increase in dark energy parameter constraints relative to the constraints 
provided by a survey with our standard tomography. In the bottom row, the value of the minimum photometric redshift of the survey runs along the horizontal 
axis while the different lines show different choices of the maximum photometric redshift as indicated. 



tive strategy for identifying and mitigating the influence of 
localized catastrophic photo-z errors. 
In § 3.4 we studied a simple strategy to limit the systemat- 



ics induced by catastrophic photo-z outliers. First, we showed 
that the statistical leverage of the highest redshift (z > 2.4) and 
lowest redshift (z < 0.3) source galaxies on dark energy con- 
straints is minimal. Eliminating all such galaxies from con- 
sideration in inferring dark energy parameters results in only a 
small increase in the statistical errors of dark energy equation 
of state constraints, but may eliminate some of the most severe 
systematic errors induced by localized catastrophic photo-z 
outliers. This implies that well-designed cuts on z ph will likely 
be a powerful and general means to mitigate systematics as- 
sociated with photo-z determination at a relatively small cost 
in statistical error. 

The published work that is most closely related to the 
present work is Bernstein & Huterer (2010 1. Our work is an 



extension and generalization of their study. Overall, we reach 
the same broad conclusions where the two studies are com- 
mensurable. In particular, we find that catastrophic errors of 
the localized variety must be controlled such that the rate of 
errors per galaxy is F cat < 10~ 3 if they are to induce tolerable 
systematic errors on dark energy parameters. 



Our work d iffers from and complements [Bernstein &| 



Huterer (20101 in several important ways. First, we have re- 



laxed the assumption that the true redshift distribution of the 
outlier population perfectly traces that of the core population 
within individual Source photometric redshift bins (see Eq. 
4]). Our treatment of photometric redshift errors is indepen- 
dent of the photometric redshift binning (as such errors would 
be in practice), while the approach of Ber nstein & Huterer] 
(20101 is limited to cases in which photometric errors both 



trace the galaxy distributions within the Source Bin and span 
the redshift range of the Source Bin. While contamination of 
the Target redshift bin is typically the larger source of induced 
systematic error, our generalization illustrates that the effects 
of modifications to the Source Bin are non-negligible and in 
some cases these offsets contribute significantly to the sys- 
tematic errors on dark energy parameters. Second, we have 
studied catastrophic errors in cases where the core photomet- 
ric redshift distribution is not perfectly calibrated. Accounting 
for uncertainty in the core distribution turns out to be quite im- 
portant: for a fixed catastrophe the magnitude of the induced 
systematic errors can vary by several orders of magnitude over 
a reasonable range of priors on the core distribution. Third, 
we have explored cases of correlated shifts in photo-z errors 
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that span multiple tomographic redshift bins (which will oc- 
cur in practice), the extreme example being the uniform error. 

We conclude our discussion section by referring to interest- 
ing, tangential results given in the appendix. In the appendix, 
we discuss the effect of different models of the nonlinear 
evolution of cosmological density perturbations on photomet- 
ric redshift calibration requirements. Weak lensing measure- 
ments take significant advantage of measurements on nonlin- 
ear scales in order to constrain cosmology. Previous work on 
the calibration of photometric redshifts has utilized the Pea- 
cock & Do dds ( 1996) formula (e.g. |Ma et al.|[2006l |Ma & 
Bernstein 2008); however, we find that using the more re- 
cent and more accurate fit of Smith et al. ( 2003| > significantly 
reduces the need for independent calibration of photometric 
redshifts. We have used the Smith et aL] ( |2003| l formula in the 
main body of this paper. We refer the reader to the Appen- 
dices for further details. 

5. SUMMARY 

We have adopted a simple, agnostic approach to estimate 
the levels at which uncalibrated photometric redshift outliers 
must be controlled to maximize the dark energy constraints 
from the weak lensing components of forthcoming imaging 
surveys such as DES, LSST, EUCLID, and JDEM. We present 
results for three fiducial imaging surveys: a relatively near- 
term DES -like survey; a future survey with a high surface 
density of galaxies but a relatively small fractional sky cover- 
age (DEEP); and a future survey with half-sky coverage and 
a lower galaxy surface density (WIDE). We considered two 
extreme cases of large, uncalibrated errors. In the case of 
a uniform photo-z catastrophe, we considered galaxies erro- 
neously assigned photometric redshifts that are unrelated to 
their true redshifts. In the case of a localized photo-z catastro- 
phe, we considered the erroneous placement of a small frac- 
tion of galaxies in some range of true redshifts at significantly 
different photometric redshifts. To be specific, we assigned 
galaxies in some range of true redshifts of width Az centered 
on a true redshift z cat to photometric redshifts near z[? at that 
differ significantly from z cat . For each type of error and sur- 
vey, we assessed the severity of the systematic errors on dark 
energy parameters that would be induced by catastrophic pho- 
tometric redshift errors. Our primary results are as follows. 

1 . A photometric redshift error of the uniform variety that 
is relevant for galaxies near the median redshift of the 
imaging survey, must be limited to a fraction of galaxies 
^cat <5x 10" 4 for DES or DEEP and F cat < 2 x 10" 4 
for WIDE, in order to induce systematic errors that are 
small compared to the statistical errors on Wq and w a . 

2. Localized catastrophic errors are most severe when they 
take some fraction of galaxies with true redshifts near 
the median survey redshift and assign them signifi- 
cantly higher or lower photo-z's. For DES, assignments 
to higher photo-z's are more severe than assignments to 
lower photo-z's while the opposite is true for WIDE and 
DEEP. However, the systematic errors induced by these 
two extremes differ by less than a factor of two in all 
cases. 

3. Limiting the fraction of galaxies exhibiting localized 
catastrophes at all redshifts to F cat < 3 x 10~ 3 for DES or 
-feat 10~ 3 f° r WIDE or DEEP will render them unim- 
portant. For localized catastrophes that occur over a 



range of true redshifts of width Az ca t near the median 
survey redshift, the fractional error rate must be con- 
trolled such thatF cat (Az cat /0.1) < 1 -3 x 10" 3 . 

4. Imperfect knowledge of the photo-z distribution for the 
core sample of galaxies loosens these requirements for 
uncalibrated catastrophic outlier control as depicted in 
Fig. [5] Roughly speaking, core calibration with spectro- 
scopic samples smaller than the statistical equivalent of 
Nspec < 10 5 leads to significantly reduced catastrophic 
error control requirements. Of course, in practice catas- 
trophic error control and core calibration will both im- 
prove as N spec increases. 

5. The statistical leverage of the highest redshift (z > 2.4) 
and lowest redshift (z < 0.3) source galaxies on dark 
energy constraints is small. Eliminating all such galax- 
ies from consideration in inferring dark energy param- 
eters results in a < 20% increase in the statistical errors 
on dark energy, but may eliminate the most severe sys- 
tematic errors induced by localized catastrophic photo- 
z outliers. 

6. In the appendix, we show that dark energy parameter 
forecasts that include photometric redshift uncertainty 
vary significantly depending upon the treatment of the 
nonlinearity in the matter power spectrum. In partic- 
ular, using the Sm ith et aT] ([2003 ) fitting form (as we 
do in the main text) leads to weaker photo-z calibration 

1996) 



requirements than does the IPeacock & Dodds 



formula upon which the results of |Ma et al.| ( |2006[ ) are 
based. The |Smith et al.| ( [2003|l formula has b een shown 
to be more accurate than Peacock & Dodds ( 1996 ) sug- 



gesting that degradation due to photo-z uncertainty may 
be less than M a et al.| ( [20 06) forecast. Only a rigorous 
numerical study can determine this definitively. 

This level of photometric redshift outlier control is chal- 
lenging in comparison to the yields of contemporary methods 
and data. Existing spectroscopic samples are not represen- 
tative of the galaxy populations that will be utilized to con- 
strain dark energy in forthcoming imaging surveys. For ex- 
ample, the Deep Extragalactic Evolutionary Probe (DEEP2) 
has a 7 0% success rate for obtaining spectroscopic redshifts 
( Cooper et al. 2006| l, where star-forming galaxies with z > 1 .4 



constitute roughly half of the failed targets (Freeman et al. 



2009). Moreover, spurious photo-z outliers persist even with 
techniques developed in conjunction with spectroscopic data 
that span the region of parameter space occupied by the pho- 
tometric sample. As a nearly contemporaneous example, 
|Nishizawa et al.| ( f2010| ) construct galaxies using simple spec- 
tral templates assuming a number of particular stellar popu- 
lations based on the COSMOS galaxy catalog. Applying Le- 
Phare 11 to their mock spectra gives offset islands in z-z ph 
space containing more than 5% of the probability, which re- 
mains true even after refining their redshift estimator. In the 
analysis of the COSMOS data, |Ilbert et al.| (|2"009 ) achieve an 
outlier rate of 0.7% for a subsample of their brightest objects 
(17.5 < /ab < 22.5). However, their outlier rate dramatically 
increases to 15.3% when they apply their photo-z techniques 
to a subsample of faint objects (22.5 < /ab < 24). An alter- 
native approach to photo-z calibration is adopted in |Cunha| 



et al. ( 2009 ), who applied their weighted training set method 



1 ' http://www.cfht.hawaii.edu/ amouts/LEPHARE/cfht_lephare/lephare.html 
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to both simulated and actual SDSS data. Their methods sub- 
stantially improve upon the ability to directly reconstruct the 
redshift distribution of a photometric sample, but errors in the 
reconstructed N(z) remain at the percent level. Thus, while 
contemporary photo-z codes do provide useful guidance, out- 
lier fractions greater than ~ 10~ 3 persist and can affect the 
dark energy program. The ability to either limit, or under- 
stand, such outlier populations significantly better than the 
current state-of-the-art will be necessary to exploit fully the 
promise of cosmic shear tomography. 
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APPENDIX 

THE NONLINEAR POWER SPECTRUM AND PHOTOMETRIC REDSHIFT CALIBRATION REQUIREMENTS 

Much of the constraining power of weak lensing surveys arises from measurements on scales where the structures causing 
the deflections are undergoing nonlinear gravitational evol ution (e.g. [ Huterer & Takada 2005 ). Restricting c onsideration to large 
scales significantly degrades cosmological constraints (e.g. |Huterer|2002[ |Huterer & Takada|2005 Zentner et al.|2008| |Schmidt| 
|2008 1 [Hearin & Zentner|2 009), so it is necessary to model nonlinear evolution in order to utilize weak lensing to constrain dark 
energy. At least three approximate and related techniques are in common use: (1) the fitting formula of Peacock & Dodds ( 1996), 
which is based on the HKLM method (|Hamilton et al.|1991); ( 2) the halo model ( |Scherrer & Bertschinger|1991[ |S~eljak 2000; 



Ma & Fry|2000l|Scoccimarro et al.|200T|TCoo ray & Sheth 2002); and (3) the fitting formula of |Smith et al.| ( |2003| l. T he works of 
al. C" ~ ' " ■ ' . . - 



Ma et al.|(|2006[ ) and |Ma & Bernste in (2008 ) specifying requirements for photometric redshift calibration employ the Peacock & 
Dodds| ( | 1996j relation. 

In the course of our study, we have recomputed the photometric redshift calibration requirements using each of the three 
approximate techniques mentioned in the previous paragraph. In the limit of perfect knowledge of the photometric redshift 
distribution, each of these fitting formulas gives nearly identical dark energy constraints. However, we have found that the 
photometric redshift calibration requirements have a strong dependence upon the method used to model nonlinear structure. We 
summarize this finding in Figure Al where we display contours of constant degradation in the statistical error on w a as a function 
of both the prior on the bias Az and the prior on the dispersion a z . In other words, we show contours of cr(w a ) in units of 
the statistical constraint on w a in the limit that the photo-z distribution parameters are known perfectly prior to the weak lensing 
analysis, <7 per f(w a ). We assume that the same priors are applied to all of our 31 dispersion parameters and 31 bias parameters at 
each redshift bin (see § 2.1.1 1. We summarize our findings in this way so that these results can be compared directly to Figure 7 
in |Ma et aT7| ( |2006) l. To make the comparison as direct as possible, we have computed these forecasts using the fiducial cosmology 
and experimental setup of |Ma et aL| ( |2006| l, which differs slightly from those considered in the main text. In this appendix only, 
our fiducial cosmology is o>m = 0.14, u>b = 0.024, n s = 1.0, A 2 ^ = 2.4 x 10~ 5 (giving a% ~ 0.91), and J7de = 0.73 combined with 

experimental parameters of / s k y = 0.1 and N A = 55 arcmin -2 . 

The upper left panel in Figure Al shows photo-z calibration requirements estimated using the Peacock & Dodds ( 1996) treat- 
ment of nonlinear power. This panel shows nearly identical results to those in Figure 7 of (Ma et al. 2006) so that this panel 
validates our methods and provides a useful baseline to compare with the other panels. According to this result, ensuring that 
constraints on w a are not degraded by more than a factor of two requires knowing the photo-z dispersion and bias to roughly 
~ 1% prior to undertaking the weak lensing analysis. The upper, right panel of Fig. |Al| shows the same requirements constructed 
using the halo model for nonlinear clustering. In the limit of restrictive prior knowledge of the photo-z distribution the Peacock 
|& Dodds] ( |1996| l and halo model results yield nearly the same constraints. When the photo-z distributions have relatively unre- 
stnctive priors, the two techniques yield moderately different levels of projected degradation with, for example, uncertainty in 
the photo-z dispersion of A(<r z ) w 1 corresponding to a factor of ten degradation in the Peacock & Dodds ( 1996) case but a factor 
of six degradation in the halo model calculation. 

The largest differences among the forecasts comes from comparing the requirements using Peacock & Dodds| ( |1996[ ) to those 
computed using the Smi th et al.| ([2003 ) fit. As with the halo model comparison, the different techniques agree well when prior 
knowledge of the photo-z distribution is very restrictive; thus as long as degradations due to photo-z uncertainty are < 10-20% 
then it does not matter which technique one uses to predict the nonlinear evolution. It is interesting that the constraints in the 
case of the halo model treatment degrade significantly less rapidly as prior knowledge becomes less and less restrictive. Turning 
to the Smith et al. ( 2003| > fit, one would conclude that ensuring less than a factor of two degradation on the w a constraint requires 
~ 18% knowledge of the dispersion and ~ 40% knowledge of the bias as compared to the ~ 1% requirements that result from 
the Peacock & Dodds| ([T996 ) analysis. 

Clearly, at most one of these treatments can represent the growth of cosmic structure faithfully. In the main text, we presented 
results using the Smith et aL] ( 2003| l formula because these authors perform a detailed numerical study that finds the Peacock & 
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FIG. Al . — Contour plots for the level of w a constraint degradation as a function of priors on the photometric redshift scatter <r z and bias z bms - In this case, the 
priors are applied uniformly to the photometric redshift parameters are each redshift. The contours demarcate equal parameter degradation defined as the error 
on vv a after marginalizing over photometric reshift uncertainties. We show constraints in units of the equivalant constraint in the limit of perfect knowledge of 
photometric redshift parameters, S = o"(w a )/o"perf(Wa). The upper, left panel was computed using the Peacock & Dodds ( 1996 1 fitting formula for the nonlinear 
power spectrum of density fluctuation s and amounts to a near reproduction of the right panel of Figure / in Ma et al. (zuuoi. The upper, right panel was computed 
using the halo model as described in Zentner et al. 1 2008 1. The bottom panel was computed using the Smith et al. ( 2003 1 relation for the nonlinear power spectrum 
of density fluctuations. Significant diff erences between the levels of degradation are evident. Note that in this figure, we use a different set of cosmological and 
experimental parameters so that this result is directly comparable to those in Figure 7 of Ma et al. 1 2006 1. 



Dodds ( 1996 1 and simple implementations of the halo model to be imprecise on scales relevant for cosmic shear cosmology. In 
the context of these fitting formulae, we find that Smith et al. ( 2003| l predicts greater power than |Peacock & Dodds ( 1996! on 
scales most relevant to lensing (0.1 < fc/ZzMpc -1 < 10), particularly at high redshift. At this point, it is not possible to make a firm 



statement as to which approach is correct, but an exhaustive simulation program similar to that being carried out by Heitmann 



et al.| ( |2005| |20 08 2009) may be capable of providing a more definitive resolution in the case of dissipationless evolution. 
Additional effort will be needed to treat any modifications induced by the baryonic component of the universe ( |White|2004| |Zhari] 
& Knox|2004"l|Jing et al.|2006l|Rudd et al.|2007||Zentner et al.|2008l|Stanek et al.|2009[[GuiiTet et al.|2009[ ). 



REFERENCES 



Abdalla, F. B., Banerji, M., Lahav, O., & Rashkov, V. 2008, ArXiv:0812.3831 
Abell, P. etal. 2009, [arXiv:0912.0201] 

Albrecht, A., Bernstein, G., Cahn, R., Freedman, W. L., Hewitt, J., Hu, W., 
Huth, J., Kamionkowski, M., Kolb, E. W., Knox, L., Mather, J. C, Staggs, 
S., & Suntzeff, N. B. 2006, (astro-ph/0609591) 

Banerji, M., Abdalla, F. B., Lahav, O., & Lin, H. 2008, MNRAS, 386, 1219 



Benjamin, J., Heymans, C, Semboloni, E., van Waerbeke, L., Hoekstra, H., 
Erben, T., Gladders, M. D., Hetterscheidt, M., Mellier, Y., & Yee, H. K. C. 
2007, MNRAS, 381,702 
Bernstein, G. & Huterer, D. 2010, MNRAS, 401, 1399 
Bolzonella, M., Miralles, J., & Pello, R. 2000, A&A, 363, 476 
Brammer, G. B., van Dokkum, P. G., & Coppi, P. 2008, ApJ, 686, 1503 



18 



Hearin et al. 



Chevallier, M. & Polarski, D. 2001, International Journal of Modern Physics 

D, 10, 213 

Collister, A. A. & Lahav, O. 2004, PASP, 116, 345 

Cooper, M. C, Newman, J. A., Croton, D. J., Weiner, B. J., Willmer, C. N. A., 
Gerke, B. R, Madgwick, D. S., Faber, S. M., Davis, M., Coil, A. L., 
Finkbeiner, D. P., Guhathakurta, P., & Koo, D. C. 2006, MNRAS, 370, 
198 

Cooray, A. & Hu, W. 2001, ApJ, 554, 56 
Cooray, A. & Sheth, R. 2002, Phys. Rep., 372, 1 

Coupon, J., Ilbert, O., Kilbinger, M., McCracken, H. J., Mellier, Y., Arnouts, 
S., Bertin, E., Hudelot, P., Schultheis, M., Le Fevre, O., Le Brun, V., 
Guzzo, L., Bardelli, S., Zucca, E., Bolzonella, M., Garilli, B., Zamorani, 

G. , Zanichelli, A., Tresse, L., & Aussel, H. 2009, A&A, 500, 981 
Cunha, C. E., Lima, M., Oyaizu, H., Frieman, J., & Lin, H. 2009, MNRAS, 

396, 2379 

Dahlen, T., Mobasher, B., Jouvel, S., Kneib, J.-P, Ilbert, O., Arnouts, S., 

Bernstein, G., & Rhodes, J. 2008, AJ, 136, 1361 
Dodelson, S., Shapiro, C, & White, M. 2006, Phys. Rev. D, 73, 023009 
Dodelson, S. & Zhang, P. 2005, Phys. Rev. D, 72, 083001 
Dore, O., Martig, M., Mellier, Y., Kilbinger, M., Benjamin, J., Fu, 

L., Hoekstra, H., Schultheis, M., Semboloni, E., & Tereno, I. 2007, 

ArXiv:0712.1599 

Feldmann, R., Carollo, C. M., Porciani, C, Lilly, S. J., Capak, P., Taniguchi, 
Y, Le Fevre, O., Renzini, A., Scoville, N., Ajiki, M., Aussel, H., Contini, 
T., McCracken, H., Mobasher, B., Murayama, T., Sanders, D., Sasaki, 
S., Scarlata, C, Scodeggio, M., Shioya, Y, Silverman, J., Takahashi, M., 
Thompson, D., & Zamorani, G. 2006, MNRAS, 372, 565 

Freeman, P. E., Newman, J. A., Lee, A. B., Richards, J. W., & Schafer, C. M. 
2009, MNRAS, 398, 2012 

Fu, L., Semboloni, E., Hoekstra, H., Kilbinger, M., van Waerbeke, L., Tereno, 
I., Mellier, Y, Heymans, C, Coupon, J., Benabed, K., Benjamin, J., Bertin, 

E. , Dore, 0., Hudson, M. J., Ilbert, O., Maoli, R., Marmo, C, McCracken, 

H. J., & Menard, B. 2008, A&A, 479, 9 

Guillet, T, Teyssier, R., & Colombi, S. 2009, ArXiv:0905.2615 
Hamilton, A. J. S., Kumar, P., Lu, E., & Matthews, A. 1991, ApJ, 374, LI 
Hearin, A. P. & Zentner, A. R. 2009, Journal of Cosmology and Astro-Particle 

Physics, 4, 32 
Heavens, A. 2003, MNRAS, 343, 1327 

Heitmann, K, Higdon, D., White, M., Habib, S., Williams, B. J., & Wagner, 

C. 2009, ArXiv:0902.0429 

Heitmann, K., Ricker, P. M., Warren, M. S., & Habib, S. 2005, ApJS, 160, 28 
Heitmann, K., White, M., Wagner, C, Habib, S., & Higdon, D. 2008, 
ArXiv:0812.1052 

Hoekstra, H. & Jain, B. 2008, Annual Review of Nuclear and Particle 
Science, 58, 99 

Hoekstra, H., Yee, H. K. C, & Gladders, M. D. 2002, ApJ, 577, 595 

Hu, W. 1999, ApJ, 522, L21 

Hu, W. & Tegmark, M. 1999, ApJ, 514, L65 

Huterer, D. 2002, Phys. Rev. D, 65, 063001 

Huterer, D. & Takada, M. 2005, Astroparticle Physics, 23, 369 

Huterer, D., Takada, M., Bernstein, G., & Jain, B. 2006, MNRAS, 366, 101 

Huterer, D. & Turner, M. S. 2001, Phys. Rev. D, 64, 123527 

Ilbert, O., Capak, P., Salvato, M., Aussel, H., McCracken, H. J., Sanders, 

D. B., Scoville, N, Kartaltepe, J., Arnouts, S., Floc'h, E. L., Mobasher, 
B., Taniguchi, Y, Lamareille, F, Leauthaud, A., Sasaki, S., Thompson, D., 
Zamojski, M., Zamorani, G., Bardelli, S., Bolzonella, M., Bongiorno, A., 
Brusa, M., Caputi, K. I., Carollo, C. M., Contini, T., Cook, R., Coppa, G., 
Cucciati, O., de la Torre, S., de Ravel, L., Franzetti, P., Garilli, B., Hasinger, 
G, Iovino, A., Kampczyk, P., Kneib, J., Knobel, C, Kovac, K., LeBorgne, 
J. E, LeBran, V., Fevre, O. L., Lilly, S., Looper, D., Maier, C, Mainieri, 
V., Mellier, Y, Mignoli, M., Murayama, T, Pello, R., Peng, Y, Perez- 
Montero, E., Renzini, A., Ricciardelli, E., Schiminovich, D., Scodeggio, 
M., Shioya, Y, Silverman, J., Surace, J., Tanaka, M., Tasca, L., Tresse, L., 
Vergani, D., & Zucca, E. 2009, ApJ, 690, 1236 

Ishak, M. 2005, MNRAS, 363, 469 

Jarvis, M., Bernstein, G. M., Fischer, P., Smith, D., Jain, B., Tyson, J. A., & 

Wittman, D. 2003, AJ, 125, 1014 
Jarvis, M., Jain, B., Bernstein, G., & Dolney, D. 2006, ApJ, 644, 71 
Jing, Y. P., Zhang, P., Lin, W. P., Gao, L., & Springel, V. 2006, ApJ, 640, 

L119 

Jouvel, S., Kneib, J. P., Ilbert, O., Bernstein, G., Arnouts, S., Dahlen, T., 
Ealet, A., Milliard, B., Aussel, H., Capak, P., Koekemoer, A., Le Brun, 
V, McCracken, H., Salvato, M., & Scoville, N. 2009, ApJ, 504, 359 



Jungman, G, Kamionkowski, M., Kosowsky, A., & Spergel, D. N. 1996, 

Phys. Rev. D, 54, 1332 
Kitching, T. D., Heavens, A. F, Taylor, A. N., Brown, M. L., Meisenheimer, 

K, Wolf, C, Gray, M. E., & Bacon, D. J. 2007, MNRAS, 376, 771 
Kitching, T. D., Taylor, A. N, & Heavens, A. F. 2008, MNRAS, 389, 173 
Komatsu, E., Dunkley, J., Nolta, M. R., Bennett, C. L., Gold, B., Hinshaw, G., 

Jarosik, N., Larson, D., Limon, M., Page, L., Spergel, D. N., Halpern, M., 

Hill, R. S., Kogut, A., Meyer, S. S., Tucker, G. S., Weiland, J. L., Wollack, 

E., & Wright, E. L. 2008, ArXiv:0803.()547 
Kosowsky, A., Milosavljevic, M., & Jimenez, R. 2002, Phys. Rev. D, 66, 

063007 

Lima, M., Cunha, C. E., Oyaizu, H., Frieman, J., Lin, H., & Sheldon, E. S. 

2008, MNRAS, 390, 118 
Lima, M. & Hu, W. 2007, Phys. Rev. D, 76, 123013 
Linder, E. V. 2003, Physical Review Letters, 90, 091301 
Ma, C.-P. & Fry, J. N. 2000, ApJ, 543, 503 
Ma, Z. & Bernstein, G. 2008, ApJ, 682, 39 
Ma, Z., Hu, W., & Huterer, D. 2006, ApJ, 636, 21 
Margoniner, V. E. & Wittman, D. M. 2008, ApJ, 679, 31 
Munshi, D., Valageas, P., van Waerbeke, L., & Heavens, A. 2008, Phys. Rep., 

462, 67 

Newman, J. A. 2008, ApJ, 684, 88 
Newman, J. A. et al. 2010, In Prepartion 

Nishizawa, A. J., Takada, M., Hamana, T., & Furusawa, H. 2010, 
ArXiv: 1002.2476 

Oyaizu, H., Lima, M., Cunha, C. E., Lin, H., & Frieman, J. 2008, ApJ, 689, 
709 

Peacock, J. A. & Dodds, S. J. 1996, MNRAS, 280, L19 

Pen, U.-L., Zhang, T., van Waerbeke, L., Mellier, Y, Zhang, P., & Dubinski, 

J. 2003, ApJ, 592, 664 
Refregier, A. 2003, ARA&A, 41, 645 

Refregier, A., Amara, A., Kitching, T. D., Rassat, A., Scaramella, R., Weller, 
J., & Euclid Imaging Consortium, f. t. 2010, ArXiv:1001.0061 

Refregier, A., Massey, R., Rhodes, J., Ellis, R., Albert, J., Bacon, D., 
Bernstein, G., McKay, T, & Perlmutter, S. 2004, AJ, 127, 3102 

Rudd, D. H., Zentner, A. R., & Kravtsov, A. V. 2007, ApJ, 672, 19 

Scherrer, R. J. & Bertschinger, E. 1991, ApJ, 381, 349 

Schmidt, F. 2008, Phys. Rev. D, 78, 043002 

Schulz, A. E. 2009, ArXiv:0910.3683 

Scoccimarro, R., Sheth, R. K., Hui, L., & Jain, B. 2001, ApJ, 546, 20 
Seljak, U. 1997, ApJ, 482, 6 
— . 2000, MNRAS, 318, 203 

Semboloni, E., Mellier, Y, van Waerbeke, L., Hoekstra, H., Tereno, I., 

Benabed, K., Gwyn, S. D. J., Fu, L., Hudson, M. J., Maoli, R., & Parker, 

L. C. 2006, A&A, 452,51 
Smith, R. E., Peacock, J. A., Jenkins, A., White, S. D. M., Frenk, C. S., 

Pearce, F. R., Thomas, P. A., Efstathiou, G., & Couchman, H. M. P. 2003, 

MNRAS, 341, 1311 
Song, Y.-S. & Knox, L. 2004, Phys. Rev. D, 70, 063510 
Stanek, R., Rudd, D., & Evrard, A. E. 2009, MNRAS, 394, Lll 
Sun, L., Fan, Z., Tao, C, Kneib, J., Jouvel, S., & Tilquin, A. 2009, ApJ, 699, 

958 

Takada, M. & Jain, B. 2004, MNRAS, 348, 897 

Takada, M. & White, M. 2004, ApJ, 601, LI 

Tegmark, M., Taylor, A. N, & Heavens, A. F. 1997, ApJ, 480, 22 

Vale, C. & White, M. 2003, ApJ, 592, 699 

Van Waerbeke, L., Mellier, Y, & Hoekstra, H. 2005, A&A, 429, 75 
White, M. 2004, Astroparticle Physics, 22, 21 1 
White, M. & Hu, W. 2000, ApJ, 537, 1 
Zentner, A. R. & Bhattacharya, S. 2009, ApJ, 693, 1543 
Zentner, A. R., Rudd, D. H., & Hu, W. 2008, Phys. Rev. D, 77, 043507 
Zhan, H. 2006, Journal of Cosmology and Astro-Particle Physics, 8, 8 
Zhan, H. & Knox, L. 2004, ApJ, 616, L75 
Zhang, P., Pen, U., & Bernstein, G. 2009, ArXiv:0910.4181 
Zhao, G.-B., Pogosian, L., Silvestri, A., & Zylberberg, J. 2009, 
ArXiv:0905.1326 



